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ABSTRACT 


In  developing  an  energy-based  design  approach  and  assessing  the  damage  potential  of 
structures,  one  must  know  the  distribution  of  earthquake  input  energy  among  energy 
components:  kinetic,  elastic  strain,  hysteretic,  and  damping.  This  report  examines  the 
influences  of  the  ground  motion  characteristics:  intensity,  frequency  content,  and  duration 
of  strong  motion  and  the  structural  properties:  ductility,  damping,  and  hysteretic  behavior 
on  the  distribution  of  input  energy  for  a one-  and  a five-story  building  using  20 
accelerograms,  ten  with  short  and  ten  with  long  duration  of  strong  motion. 

Results  indicate  that  for  certain  damping  ratios,  ductility  has  a significant  influence  on 
input  energy  and  its  distribution  among  energy  components  in  a structure.  For  a given 
ductility  ratio,  small  damping  ratio  (less  than  5%)  has  a minor  effect  on  input  energy,  but 
a major  influence  on  the  energy  distribution.  Damping  ratios  larger  than  5%  have  a 
significant  influence  on  the  input  energy  and  its  distribution.  Three  energy  ratios  that 
relate  to  hysteretic  energy  were  computed:  the  maximum  ratio  of  hysteretic  to  input 
energy  (Eh/Eir)m,  the  ratio  of  the  maximum  hysteretic  energy  to  the  maximum  input 
energy  Ehm/Einn , and  the  equivalent  number  of  yield  excursions  Neq=EhJ(Fy.up)  where  Fy 
is  the  yield  strength,  and  up  is  the  plastic  deformation.  It  is  found  that  ( Ef,/Eir)„ , generally 
reflects  the  energy  demand  for  the  largest  yield  excursion,  and  EhJEirm  and  Neq  reflect  the 
energy  demand  for  the  entire  duration  of  accelerogram.  The  study  shows  that  ( Eh/E,r)m  is 
independent  of  the  duration  of  strong  motion  and  period  of  structure;  however,  EhmIEirm  is 
independent  of  both  only  for  periods  less  than  1 s.  Results  indicate  that  as  the  duration 
becomes  longer  the  equivalent  number  of  yield  excursions  Neq  increases  indicating  more 
structural  damage. 

The  influence  of  ground  motion  characteristics  and  structural  properties  on  the 
distribution  of  energy  parameters  for  a five-story  building  with  fixed-base,  base-isolation, 
supplemental  damping,  and  semi-active  control  are  examined  using  the  20  accelerograms. 
The  results  show  that:  1)  the  distribution  of  energy  through  the  height  of  the  building  is 
mostly  independent  of  the  frequency  content  and  the  duration  of  strong  motion,  2)  base- 
isolation,  supplemental  damping,  and  semi-active  control  reduce  the  damage  potential  by 
reducing  the  input  and  hysteretic  energy  demands  and  have  significant  influences  on  the 
distribution  of  energy  through  the  height  of  the  building. 
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1. 


INTRODUCTION 


The  current  seismic  design  practice  which  is  based  on  strength  principles  (using  the 
acceleration  spectra)  does  not  directly  account  for  the  influence  of  the  duration  of  strong 
motion  nor  for  the  hysteretic  behavior  of  the  structure.  The  hysteretic  behavior  is 
addressed  indirectly  by  using  the  response  modification  factor  R which  is  based  primarily 
on  the  structural  system  selected.  A design  approach  based  on  energy  (Vision  2000),  on 
the  other  hand,  has  the  potential  to  address  the  effects  of  the  duration  and  hysteretic 
behavior  directly. 

The  earthquake  input  energy  transmitted  to  a structure  consists  of  the  kinetic  energy, 
elastic  strain  energy,  damping  energy,  and  hysteretic  energy.  Kinetic  energy  reflects  the 
work  of  the  inertia  force.  Elastic  strain  energy  is  the  portion  of  the  input  energy  stored  in 
the  structure  in  the  form  of  elastic  strain.  Damping  energy  is  the  work  of  the  damping 
force.  Hysteretic  energy  is  the  energy  dissipated  through  the  hysteretic  action  and  is 
associated  with  the  damage  potential  of  the  structure  (Kuwamura  and  Galambos,  1989). 

Input  energy  reflects  the  intensity  of  earthquake  ground  motion.  It  is  a more  appropiate 
measure  of  the  intensity  than  the  peak  ground  acceleration  or  the  effective  peak 
acceleration  defined  by  ATC  3-06  (1978)  as 


EPA  = SJ  2.5 


(1.1) 


where  Sa  is  the  average  spectral  acceleration  in  the  constant  acceleration  region,  or  the 
Spectrum  Intensity  defined  by  Housner  (1975)  as 


SI=  \~PSVdT 
Jo.i  " 


(1.2) 


where  P SV  is  the  pseudo  spectral  velocity  and  T„  denotes  the  period  of  structure  in 
seconds.  Earthquake  input  energy  accounts  for  both  the  ground  motion  characteristics 
(frequency  content,  intensity,  and  duration  of  strong  motion)  and  structural  properties 
(ductility,  damping,  period,  and  hysteretic  behavior),  whereas  PGA , EPA.  and  SI  account 
for  some  but  not  all  these  characteristics. 

A design  approach  based  on  energy  reflects  the  distribution  of  input  energy  and  is  suited 
for  seismic  design  of  structures  with  base  isolation  system  and/or  supplemental  damping 
devices.  In  an  energy-based  design  approach,  once  the  energy  demand  for  a structure  is 
estimated  from  the  earthquake  ground  motion,  the  damage  potential  can  be  quantified  by 
a combination  of  response  and  energy  parameters  (see  Park  and  Ang,  1985).  Sufficient 
strength  and  energy  dissipation  capacity  should  be  provided  in  the  structure  for  an 
acceptable  damage  threshold,  i.e.  a desired  performance  objective. 
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In  seismic  design  of  structures,  the  strength  demand  is  defined  by  the  story  shear  and  the 
energy  demand  should  be  defined  by  the  hysteretic  energy  for  each  story.  While  the 
strength  and  energy  dissipation  capacities  can  be  increased  by  using  members  with  larger 
area  and  section  modulus  or  by  using  materials  with  greater  yield  strength,  the  energy 
dissipation  capacity  can  be  increased  by  using  ductile  structural  systems  such  as  the 
eccentrically  braced  frames  (EBF)  instead  of  concentrically  braced  frames  (CBF),  see 
Newmark  and  Rosenbluth  1971.  The  energy  dissipation  capacity  for  moment  resisting 
frames  can  be  increased  by  providing  special  detailing  in  steel  and  reinforced  concrete 
buildings  (refer  to  IBC  2000).  Computing  the  energy  dissipation  capacity  of  various 
structural  members  is  not  within  the  scope  of  this  report  (one  should  refer  to  Uang  and 
Bertero,  1990).  Only  the  energy  demand  is  addressed  herein. 

This  report  examines  the  influences  of  ground  motion  characteristics  and  structural 
properties  on  earthquake  input  energy  and  its  distribution  among  energy  components  for 
a one-  and  a five-story  building  using  20  accelerograms,  ten  with  short  and  ten  with  long 
duration  of  strong  motion.  A summary  of  previous  research  on  this  subject  is  presented 
in  Chapter  2.  Chapter  3 presents  the  energy  balance  equation  and  the  energy  time- 
histories  for  an  SDOF  structure.  Chapter  4 examines  the  existing  procedures  for 
computing  the  earthquake  input  energy,  and  how  ground  motion  characteristics  and 
structural  properties  influence  the  input  energy.  In  Chapter  5,  the  relationship  between 
the  hysteretic  energy  and  damage  potential  is  examined.  The  distribution  of  story  shear 
and  energy  in  a 5-story  building  without  control  devices,  with  base  isolation  system,  with 
supplemental  damping  devices,  and  with  semi-active  control  are  presented  in  Chapters  6 
and  7.  Finally,  the  conclusions  are  discussed  in  Chapter  8. 
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2. 


SUMMARY  OF  PREVIOUS  RESEARCH 


McKevitte  et  al.  (1980)  computed  the  input  energy,  hysteretic  energy,  and  the  ratio  of  the 
cumulative  hysteretic  energy  to  input  energy  for  SDOF  and  MDOF  structures  (three-  and 
ten-story)  with  different  structural  properties  subjected  to  four  earthquake  records  (El 
Centro  1940,  Taft  1952,  Parkville  1956,  and  Pacoima  Dam  1971).  They  concluded  that 
the  energy  dissipated  through  inelastic  deformation  depends  on  the  force-deformation 
characteristics,  yield  strength,  and  damping.  They  observed  that  the  percentage  of  input 
energy  dissipated  by  the  hysteretic  action  was  approximately  the  same  for  different 
records.  McKevitte  et  al.  concluded  that  the  ratio  of  the  maximum  hysteretic  energy  to 
the  maximum  input  energy  for  an  MDOF  structure  can  be  estimated  from  an  SDOF 
structure  with  the  same  fundamental  period,  yield  strength,  and  damping. 

Zahrah  and  Hall  (1984)  computed  the  input  energy  for  eight  earthquake  records  and 
proposed  an  equivalent  number  of  yield  excursions  to  quantify  the  earthquake  damage 
potential.  They  observed  that  ductility,  damping,  and  the  post-  to  pre-yield  stiffness 
ratios  have  small  effects  on  the  input  and  hysteretic  energies  for  a structure  with  bilinear 
behavior.  They  stated  that  the  equivalent  number  of  yield  excursions  may  be  different  for 
different  accelerograms. 

Akiyama  (1985)  compared  the  input  energy  from  the  Fourier  Spectra  of  ground 
acceleration  for  a five-story  building  with  different  structural  properties,  and  for  an 
equivalent  one-story  building  having  the  same  fundamental  period,  total  mass,  and  yield 
strength  using  the  S00E  component  of  the  1940  El  Centro  record.  He  showed  that  the 
total  input  energy  transmitted  to  a five-story  building  can  be  computed  from  the  input 
energy  transmitted  to  the  equivalent  one-story  building  with  the  same  fundamental 
period,  damping,  and  yield  strength. 

Nakashima  et  al.  ( 1 996)  investigated  the  effect  of  large  post-  to  pre-yield  stiffness  ratios 
(up  to  0.75)  for  bilinear  SDOF  and  MDOF  structures  (5,  10,  13,  15,  17,  and  20  story) 
subjected  to  three  long  duration  records  (El  Centro  1940,  Taft  1952,  and  Miyagiken-oki 
1978).  They  concluded  that  in  general,  large  post-  to  pre-yield  stiffness  ratios  have  a 
minor  effect  on  the  input  energy;  however,  as  the  post-  to  pre-yield  stiffness  ratio 
increases,  the  input  energy  spectra  increase  in  the  vicinity  of  the  predominant  period. 
They  also  concluded  that  the  energy  transmitted  to  an  MDOF  structure  can  be  computed 
from  the  input  energy  transmitted  to  an  equivalent  SDOF  structure. 

Rahnama  and  Manuel  (1996)  computed  the  input  energy,  the  hysteretic  energy,  and  the 
ratio  of  the  cumulative  hysteretic  energy  to  input  energy  for  ductility  ratios  of  2 and  6 
with  5%  damping  for  six  sets  of  19  accelerograms  each:  actual  records,  simulated  records 
with  the  same  duration  as  the  actual  records,  and  four  sets  of  simulated  records  with 
durations  of  5,  10,  15,  and  20  s.  They  concluded  that  as  the  duration  increases,  the  input 
and  hysteretic  energies  increase.  The  duration,  however,  does  not  influence  the  ratio  of 
cumulative  hysteretic  energy  to  input  energy. 


3 


Goel  (1997)  investigated  the  distribution  of  energy  in  asymmetric  structures.  He 
concluded  that  the  input  energy  transmitted  to  a symmetric  or  to  an  asymmetric  one-story 
building  is  approximately  the  same.  He  also  showed  that  in  an  asymmetric  building  the 
hysteretic  energy  demand  is  larger  for  flexible  elements,  whereas  it  is  approximately  the 
same  for  stiff  elements. 
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3.  EARTHQUAKE  INPUT  ENERGY 


Uang  and  Bertero  (1990)  proposed  two  procedures  for  computing  the  earthquake  input 
energy:  one  based  on  the  absolute  motion  and  the  other  on  relative  motion.  The 
difference  between  the  two  procedures  is  less  important  in  damage  assessment  since  the 
hysteretic  energy,  which  is  associated  with  the  damage  potential  of  structures,  is 
independent  of  the  approach  used.  Chopra  (1995),  and  Bruneau  and  Wang  (1996) 
believe  that  the  input  energy  in  terms  of  the  relative  motion  is  more  meaningful  than  the 
input  energy  in  terms  of  the  absolute  motion  since  internal  forces  within  a structure  are 
computed  using  relative  displacements  and  velocities.  Therefore,  the  procedure  based  on 
the  relative  displacement  is  used  herein. 


3.1.  Energy  Balance  Equation 

The  energy  balance  equation  for  an  SDOF  structure  can  be  written  as  (Uang  and  Bertero 
1990) 


\ m ii  u dr  + \ c u2  dr  + [ fii  dr  - - [ mil  ii  dr 

Jo  Jo  Jo  5 Jo  8 


(3.1) 


where  m is  mass  of  the  structure,  c is  the  damping  coefficient,  Js  is  the  restoring  force,  u 
is  the  relative  displacement  of  the  mass  with  respect  to  ground,  ii  is  the  velocity  of  the 
mass  with  respect  to  ground,  ii  is  the  acceleration  of  the  mass  with  respect  to  ground, 
denotes  the  ground  acceleration,  and  t is  time.  The  energy  balance  equation  (3.1)  can 

be  written  as 


Eh.  + Ed+Ea=Ew  (3.2) 

where  Ej^  is  the  relative  kinetic  energy,  Ed  is  the  damping  energy,  Ea  is  the  absorbed 
energy,  and  Eir  denotes  the  relative  input  energy  such  that 


Ekr  = Vi  mu2 

(3.3) 

ej  = {,/'*' 2 dt 

(3.4) 

E«  = I dT 

(3.5) 

Eir——[  miiiidr 

Jo  8 

(3.6) 
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The  absorbed  energy  consists  of  the  recoverable  elastic  strain  energy  Es  and  the 
irrecoverable  hysteretic  energy  Eh  where 


(3.7) 

(3.8) 

in  which  k is  the  pre-yield  stiffness  of  the  structure. 

Figure  3.1  shows  the  energy  time-histories  for  a bilinear  SDOF  structure  with  a pre-yield 
period  T„  = 0.5  5,  a ductility  ratio  fj  = 3,  a pre-yield  damping  ratio  £=  5%,  and  a post-  to 
pre-yield  stiffness  ratio  a = 0.02  subjected  to  the  S00E  component  of  El  Centro,  the 
Imperial  Valley,  California  earthquake  of  May  18,  1940.  The  Bouc-Wen  model  (Bouc, 
1967;  Wen,  1976)  is  used  to  characterize  the  nonlinear  behavior  of  the  structure,  see 
Appendix  A.  Parameters  used  for  characterizing  the  bilinear  behavior  with  Bouc-Wen 
model  are:  A - 1,  /?=  y-  0.5,  n = 20,  and  a = 0.02.  To  overcome  the  problem  associated 
with  the  non-zero  initial  conditions  of  the  records,  particularly  in  non-linear  analyses  (see 
Pecknold  and  Riddell,  1978),  a one-second  acceleration  pulse  is  added  to  the  beginning 
of  the  accelerogram. 

Figure  3.1  shows  that  the  maximum  kinetic  and  elastic  strain  energies  occur  in  the  initial 
stages  of  the  excitation,  whereas  the  maximum  damping,  hysteretic,  and  input  energies 
occur  at  the  end  of  the  excitation.  Therefore,  the  duration  of  strong  motion  significantly 
affects  the  maximum  damping  energy,  the  maximum  hysteretic  energy,  and  the 
maximum  input  energies,  but  not  the  maximum  kinetic  energy  and  the  maximum  elastic 
strain  energy. 


Es=fs2/2k 
E„  = Ea  - Es 


3.2.  Energy  Ratio  Time-Histories 

Figure  3.2  shows  the  time-histories  of  the  ratio  of  energy  components  to  input  energy  for 
a bilinear  SDOF  structure  with  a pre-yield  period  T„  = 0.5  s , a ductility  ratio  /i=  3,  and  a 
5%  pre-yield  damping  ratio  subjected  to  the  S00E  component  of  the  1940  El  Centro 
accelerogram.  The  figure  indicates  that  the  time-histories  of  the  Ek,/Eir  and  Es/Eir  ratios 
are  similar  and  have  large  oscillations  during  the  linear  portion  of  the  response.  The  peak 
Ekr!E,r  and  EJEn  ratios  decay  rapidly  as  the  structure  experiences  nonlinear  deformation 
because  a significant  portion  of  the  input  energy  is  distributed  among  hysteretic  and 
damping  energies  rather  than  kinetic  and  elastic  strain  energies.  The  maximum  EhlE„- 
ratio  occurs  during  the  largest  yield  excursion  where  approximately  45%  of  the  maximum 
hysteretic  energy  Ehm  is  dissipated,  see  Figure  3.3.  The  hysteretic  energy  dissipated  in  the 
second  yield  excursion  is  32%  of  Eh,?,. 

Figures  3.2  and  3.3  indicate  that  the  maximum  ratio  of  the  hysteretic  to  input  energy 
(EJEi,)m  and  the  ratio  of  the  maximum  hysteretic  energy  to  the  maximum  input  energy 
EhJEirm  are  different,  i.e.  (EJEir)m=  11%  and  E[„JE,rm=55%.  The  (EJE^m  ratio  reflects 


6 


Figure  3.1.  Energy  time-histories  for  a bilinear  SDOF  structure  with  T„  = 0.5  s,  fi  = 3.  a = 0.02.  and 
£ = 5%  subjected  to  the  S00E  component  of  El  Centro,  the  Imperial  Valley,  CA  earthquake  of  May 
18, 1940 
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Figure  3.2.  Time-histories  of  the  ratio  of  energy  components  to  relative  input  energy  for  a bilinear 
SDOF  structure  with  T„  = 0.5  s,  a ductility  ratio  //  = 3,  and  5%  damping  subjected  to  the  S00E 
component  of  El  Centro,  the  Imperial  Valley,  CA  earthquake  of  May  18, 1940.  Parameters  for  Bouc- 
Wen  Model  are  A = 1,  f}=  y=  0.5,  n = 20,  a = 0.02 
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Figure  3.3.  Time-histories  of  Eh/Ehm  and  EhIEir  for  a bilinear  SDOF  structure  with  T„  = 0.5  s,  a 
ductility  ratio  //=  3,  and  5%  damping  subjected  to  the  S00E  component  of  El  Centro,  the  Imperial 
Valley,  CA  earthquake  of  May  18,  1940.  Parameters  for  Bouc-Wen  Model  are  A = 1,  f}=  y=  0.5,  //  = 
20,  a =0.02 
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the  damage  potential  associated  with  the  largest  yield  excursion,  and  the  EhnJE}rm  ratio 
reflects  the  damage  potential  associated  with  the  total  number  of  yield  excursions  during 
the  excitation.  The  differences  of  the  two  ratios  are  examined  in  Chapter  5. 

In  the  1994  Northridge  earthquake,  more  than  150  buildings  with  steel  moment  resisting 
frames  experienced  brittle  fracture  at  the  welded  connections,  and  in  some  cases  damage 
occurred  during  the  first  few  seconds  of  the  excitation  (SAC  95-06).  Damage 
experienced  in  the  initial  stages  of  the  excitation  may  have  been  caused  by  the 
insufficient  energy  dissipation  capacity  during  the  largest  yield  excursion  (usually  the 
first  two  yield  excursions).  Hence,  in  developing  a damage  index,  both  the  damage 
potential  for  the  largest  yield  excursion  and  for  the  entire  record  should  be  addressed. 
Although  proposing  a damage  index  is  not  within  the  scope  of  this  report  (refer  to  Park 
and  Ang,  1985),  the  relationship  between  the  hysteretic  energy  and  damage  potential  is 
discussed  in  Chapter  5. 
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4.  ESTIMATING  THE  EARTHQUAKE  INPUT  ENERGY 


In  an  energy-based  seismic  design,  one  needs  to  estimate  the  input  energy  in  a structure 
and  distribute  it  to  various  structural  components.  For  an  acceptable  damage  level,  vision 
2000  (SEAOC  1995)  and  Park  et  al.  (1987)  recommend  estimating  the  strength  and 
hysteretic  energy  demands  and  providing  sufficient  strength  and  energy  dissipation 
capacity  to  withstand  the  design  ground  motion. 

This  Chapter  presents:  1)  the  existing  methods  for  computing  the  earthquake  input 
energy,  2)  the  influence  of  ground  motion  characteristics  on  the  input  energy,  and  3)  the 
influence  of  structural  properties  on  input  energy. 

4.1  Existing  Procedures  for  Computing  Input  Energy 

Housner  (1956)  computed  the  input  energy  per  unit  mass  as 

^ = V2(PSV)2  (4.1) 

m 

where  m is  the  mass  and  PSV  denotes  the  pseudo-spectral  velocity.  He  used  equation 
(4.1)  for  elastic  and  inelastic  behavior. 

Zahrah  and  Hall  (1984)  computed  the  input  energy  per  unit  mass  as 

E r 

— - = - ii  u dr  (4.2) 

m Jo  g 

which  is  the  same  as  equation  (3.6). 

Akiyama  (1985)  proposed  the  input  energy  per  unit  mass  for  an  elastic  SDOF  structure  as 

— = '/<VEf  (4.3) 

m 

where  VE  is  an  equivalent  velocity.  He  recommended  the  following  values  for  VE  (in 
m/s): 

Ve  = 2.5T„  for  T„  < Tc 

VE  = 2.5  Tg  for  Tn>Tc  (4.4) 
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where  TG  is  the  predominant  period  of  ground  motion  as  a function  of  soil  type.  The 
values  of  TG  are  0.4,  0.6,  0.8,  and  1.0  5 for  soil  types  I (bed  rock),  II,  III,  and  IV  (softest 
soil),  respectively. 

Kuwamura  and  Galambos  (1989)  used  the  equation  proposed  by  Akiyama  (equation  4.3) 
and  recommended  the  following  values  for  Ve 


VF  = - T for  T < Tr 

2\Tg 

Ve-\^Te  for  T >Tg  (4.5) 

where  I E is  the  integral  of  the  square  of  the  ground  acceleration  for  the  total  duration  of 
accelerogram  tf 


IE  = \''usldt  (4.6) 

Using  40  accelerograms  Fajfar  et  al.  (1989)  computed  the  earthquake  input  energy  for  the 
intermediate-periods  (constant-velocity  region)  with  5%  damping  and  strength  ratios  rj  = 
0.5- 1.0,  where  r/  is  the  ratio  of  yield  force  Qy  to  mPGA  , where  PGA  is  the  peak  ground 
acceleration.  They  recommended  to  compute  the  input  energy  per  unit  mass  as 

— = 2.2((„)0!  (PGV)2  (4.7) 

m 

in  which  ta  is  the  duration  of  strong  motion  defined  by  Trifunac  and  Brady  (1975),  and 
PGV  is  the  peak  ground  velocity.  They  did  not  propose  any  formula  for  short  and  long 
periods,  i.e.  in  the  constant-acceleration  and  constant-displacement  regions. 

As  stated  previously,  Uang  and  Bertero  (1990)  proposed  two  procedures  for  computing 
the  input  energy  per  unit  mass:  one  based  on  the  absolute  motion  given  by 

^-=\u,udr  (4.8) 

m Jo  g 

where  i\t  denotes  the  total  acceleration  ( iit  = u + iig ),  and  the  other  based  on  relative 
motion  given  by 


(4.9) 
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Using  five  accelerograms  (El  Centro  1940,  Taft  1952,  Pacoima  Dam  1971,  Mexico  City 
1985,  and  San  Salvador  1986)  they  proposed  the  input  energy  for  an  SDOF  structure  with 
a ductility  ratio  ju=  5 and  a 5%  damping  ratio  be  computed  from 

— = ftl  + 0 .\2td,f(PGV)2  (4.10) 

m 2 

Kuwamura  et  al.  (1993)  computed  the  equivalent  velocity  VE  in  equation  (4.3)  as 


VE  =ij-2£iis  udt  (4.11) 

It  should  be  noted  that  the  energy  expression  under  the  square  root  sign  is  always 
positive.  Using  a different  set  of  five  accelerograms  (El  Centro  1940,  Parkfield  1966, 
Hachinohe  1968,  Pacoima  Dam  1971,  and  SCT-Mexico  1985)  they  showed  that  for  an 
undamped  elastic  structure,  the  equivalent  velocity  VE  in  equation  (4.1 1)  can  be  estimated 
by  smoothing  Vi  ( FS ): , where  FS  is  the  Fourier  spectra  given  by 


VE=FS(co)  = 


2 


Coscot  dt 


+ 


p i f 

u„  Sin  cot  dt 

Jo  8 


(4.12) 


in  which  co  is  the  circular  frequency.  Equation  (4.1 1)  used  by  Kuwamura  and  Galambos 
(1989),  and  equation  (4.2)  used  by  Zahrah  and  Hall  (1984)  are  based  on  relative  motion, 
see  equation  (4.9). 

Finally,  Manfredi  (2001)  using  244  accelerograms  examined  the  relation  among  the 
following  three  parameters: 

1)  the  equivalent  number  of  yield  excursions  N defined  by  Zahrah  and  Hall  (1984)  as 


*,= 


J hw 


J hni 


' hm 


FMm~Uy)  FyU(jH-\)  ku~(/J- 1) 


(4.13) 


where  um  and  uv  are  the  maximum  and  yield  deformations,  respectively,  and  /j.  = um  / uy 
is  the  ductility  ratio3, 

2)  the  response  modification  factor  R , and 

3)  the  dimensionless  seismic  index  Id  given  by 


( PGA)(PGV ) 


(4.14) 


a According  to  Zahrah  and  Hall  (1984)  the  minimum  value  of  Neq  is  1 .0. 
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where  Ie  is  given  by  equation  (4.6).  For  ductility  ratios  greater  than  2,  Manfredi 
proposed  to  estimate  the  input  energy  spectra  in  the  constant- velocity  region  by 


— = 0 45  ( — + 0 23  Id)(PSV)2 


(4.15) 


m 


where  //  is  the  ductility  ratio,  and  PSV  denotes  the  pseudo-spectral  velocity.  Equations 
(4.6,  4.14,  and  4.15)  show  that  as  the  duration  of  strong  motion  increases,  the  input 
energy  also  increases.  Equation  (4.15)  is  similar  to  equation  (4.1)  where  the  term  in  the 
parentheses  reflects  the  influence  of  the  duration  of  strong  motion.  Equation  (4.15)  also 
shows  that  as  the  ductility  increases,  the  input  energy  decreases. 

4.2.  The  Influence  of  Ground  Motion  Characteristics  on  Energy 


Spectra 


The  influences  of  the  earthquake  intensity,  the  duration  of  strong  motion,  and  the 
frequency  content  on  earthquake  input  energy  are  examined  using  10  accelerogrmas  with 
short  duration  (tj,  shorter  than  8 s)  and  10  with  long  duration  ( t ^ longer  than  18  s)  of 
strong  motion,  see  Table  4.1.  The  duration  td\  is  computed  using  the  definition  proposed 
by  Trifunac  and  Brady  (1975),  known  as  the  intensity-based  duration.  They  defined  the 
duration  as  the  time  interval  between  the  five  and  ninety-five  percent  contributions  to  the 
integral  of  the  square  of  ground  acceleration,  see  equation  (4.6). 

Figure  4.1  shows  the  relative  input  energy  spectra  for  a ductility  ratio  of  3 and  a 5% 
damping  ratio  for  the  S00E  component  of  El  Centro  with  effective  peak  accelerations 
EPA  = 0.3 g,  0.4g,  and  0.6g,  where  EPA  is  the  average  of  the  acceleration  spectra  in  the 
constant-acceleration  region  divided  by  2.5  (ATC  3-06,  1978),  see  equation  (1.1). 
According  to  the  figure,  as  the  effective  peak  acceleration  increases,  the  input  energy 
increases;  indicating  that  the  input  energy  is  related  to  the  intensity  of  ground  motion. 
The  energy  ratios  such  as  the  maximum  hysteretic  energy  to  the  maximum  input  energy 
EhJE,rm  (see  Figure  4.2)  are  not  affected  by  EPA;  indicating  that  the  scaling  of 
accelerograms  does  not  influence  the  distribution  of  the  earthquake  input  energy  among 
the  various  energy  components. 

Figure  4.3  presents  the  relative  input  and  hysteretic  energy  spectra  for  the  N65W 
component  of  Temblor  (tj,  = 6 s)  and  the  S69E  component  of  Taft  (tj,  — 29  s) 
accelerograms  scaled  to  an  EPA  - OAg  for  /j  = 1,3,  and  5,  and  a 5%  damping  ratio.  Both 
accelerograms  have  the  same  predominant  period  Tpe  = 0.44  s , but  the  duration  of  strong 
motion  for  the  Taft  accelerogram  is  approximately  5 times  longer  than  the  duration  for 
the  Temblor  accelerogram.  Figure  4. 3. a shows  that  ductility  has  a minor  effect  on  the 
input  and  hysteretic  energies  for  both  Taft  and  Temblor  accelerograms.  Figure  4.3 
indicates  that  on  the  average  for  //  = 1,  3,  and  5,  the  energy  spectra  for  the  Taft 
accelerogram  is  approximately  8 times  larger  than  the  energy  spectra  for  Temblor.  The 
hysteretic  energy  demand  for  the  Taft  accelerogram  is  approximately  5 and  6 times  larger 
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Table  4.1.  Strong  motion  records  used  in  the  analysis 
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Table  4.1.  Continued 
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tdi  is  the  duration  of  strong  motion  defined  by  Trifunac  and  Brady  (1975). 

Tpa  and  Tpe  are  the  predominant  periods  based  on  the  elastic  acceleration  and  input  energy  spectra,  see  Figure  4.4. 
Surface  magnitude  Ms. 


Figure  4.1.  Relative  input  energy  for  ductility  ratio  fi  = 3 and  a damping  ratio  £=  5%  for  the  SOOE 
component  of  El  Centro,  the  Imperial  Valley,  CA  earthquake  of  May  18,  1940  scaled  to  EPA  = 0.3g, 
0.4g,  and  0.6g.  Parameters  for  Bouc-Wen  Model  are^4  = 1,  f}=  y-  0.5,  n = 20,  a = 0.02 
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Figure  4.2 . Ratio  of  the  maximum  hysteretic  energy  to  the  maximum  input  energy  for  ductility  ratio 
p = 3 and  a damping  ratio  £ = 5%  for  the  SOOE  component  of  El  Centro,  the  Imperial  Valley,  CA 
earthquake  of  May  18,  1940  scaled  to  EPA  = 0.3g,  0.4g,  and  0.6g.  Parameters  for  Bouc-Wen  Model 
are  A = 1,  fi  = y = 0.5,  n = 20,  a = 0.02 
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Figure  4.3.  (a)  Earthquake  input  energy  Eir  and  (b)  hysteretic  energy  Eh  for  the  N65\N  component  of 
Temblor,  Park-field,  CA  earthquake  of  June  27,  1966  (tdi  = 6 s)  and  the  S69E  component  of  Taft 
Lincoln  School  Tunnel,  Kern  County,  CA  earthquake  of  July  21,  1952  ( tdi  = 29  s)  for  //=  1,3,  and  5, 
and  a damping  ratio  £ = 5%  scaled  to  an  EPA  = 0.4g.  Both  accelerograms  ha\e  the  same 
predominant  period  Tpe  = 0.44  s.  Parameters  for  Bouc-Wen  Model  are:  A = 1,  fi  = y = 0.5.  n = 20,  and 
a = 0.02.  The  hysteretic  energy  for  //=  1 is  not  shown  in  Figure  4.3.b  since  it  is  null. 
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than  the  hysteretic  energy  demand  for  the  Temblor  accelerogram  for  fj.  = 3 and  5, 
respectively.  Figure  4.3.b  indicates  that  as  the  duration  of  strong  motion  increases,  the 
hysteretic  energy  demand  also  increases. 

For  the  20  accelerograms  used,  the  predominant  periods  Tpe  computed  from  the  input 
energy  spectra  are  larger  than  Tpa  computed  from  the  acceleration  spectra,  see  columns  9 
and  10  of  Table  4.1.  Predominant  periods  Tpa  and  Tpe  primarily  depend  on  the  soil 
conditions.  The  average  TpeITpa  ratio  for  the  twenty  accelerograms  is  3.3.  Figure  4.4 
shows  the  normalized  acceleration  and  input  energy  spectra  for  three  accelerograms  with 
short  duration  of  strong  motion  (Figures  4.4. a - 4.4.c),  and  three  accelerograms  with  long 
duration  (Figures  4.4.d  - 4.4.ff  Since  Tpa  is  generally  different  from  Tpe,  the  strength  and 
energy  demands  do  not  necessarily  occur  at  the  same  period.  Therefore,  both  strength 
and  energy  demands  should  be  examined  in  design. 

Figure  4.5  shows  the  elastic  input  energy  spectra  for  the  90  component  of  Sylmar  and  the 
S16E  component  of  Pacoima  Dam  records  scaled  to  an  EPA  = 0.4g.  Both  accelerograms 
have  the  same  duration  of  strong  motion  (tj,  = 7 s).  The  bracketed  duration  of  the  scaled 
accelerograms  are  also  close  ( tdb  = 12  s for  Sylmar  accelerogram  and  tdb  - 11  s for 
Pacoima  Dam  accelerogram).  Figure  4.5  indicates  that  frequency  content  has  a 
significant  influence  on  input  energy  spectra,  even  for  the  same  intensity  and  duration  of 
strong  motion.  Figures  4.2  - 4.5  conclude  that  for  a given  effective  peak  acceleration  the 
influence  of  the  duration  of  strong  motion  on  the  input  energy  is  as  significant  as  the 
influence  of  the  frequency  content. 

4.3.  The  Influence  of  Structural  Properties  on  Energy  Spectra 

McKevitt  et  al.  (1980),  Zahrah  and  Hall  (1984),  Akiyama  (1985),  and  Nakashima  et  al. 
(1996)  believe  that  ductility  and  damping  do  not  have  a significant  influence  on  the 
earthquake  input  energy.  It  should  be  noted  that  these  studies  used  4,  8,  1.  and  3 
accelerograms,  respectively.  Housner  (1956),  and  Berg  and  Thomaides  (1966)  believe 
that  in  designing  a structure  to  satisfy  the  energy  demand,  the  elastic  input  energy 
computed  from  V:m(PSV)~  can  be  used  conservatively  in  lieu  of  inelastic  input  energy. 
Kato  and  Akiyama  (1982),  and  Leelataviwat  et  al.  (1999)  used  Housner’ s 
recommendation,  equation  (4.1),  to  develop  an  energy-based  design  method  for  steel 
buildings.  This  section  examines  the  influences  of  structural  properties  on  the  energy 
spectra. 

Figures  4. 6. a - 4.6.c  and  4.6.d  - 4.6.f  show  the  relative  input  energy  E,r,  equation  (4.9), 
for  three  accelerograms  with  short  and  three  accelerograms  with  long  duration  of  strong 
motion  for  ductility  ratios  ju  = 1,  2,  3,  and  5,  and  a damping  ratio  £=  5%.  It  is  observed 
that  ductility  has  a significant  influence  on  the  magnitude  and  frequency  content  of  the 
input  energy  spectra  for  the  1994  Sylmar,  1971  Pacoima  Dam,  and  1985  Mexico  City 
accelerograms  (Figures  4.6.a,  4.6.c,  and  4.6.f),  particularly  in  the  vicinity  of  the 
predominant  periods  Tpe  (definition  of  Tpe  is  illustrated  in  Figure  4.4).  Although  the 
elastic  input  energy  E„  does  not  provide  a conservative  estimate  of  the  inelastic  input 
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Figure  4.4.  Predominant  periods  Tpa  and  Tpe  from  the  elastic  acceleration  and  input  energy  spectra, 
respectively,  for  three  accelerograms  with  short  duration  of  strong  motion  (a  - c),  and  three 
accelerograms  with  long  duration  (d  - f)  in  Table  4.1.  Parameters  for  Bouc-Wen  Model  are:  A - 1, 

= y=  0.5,  n = 20,  and  a = 0.02 
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Figure  4.5.  Elastic  input  energy  spectra  for  the  90°  component  of  Sylmar  County  Hospital  Parking 
Lot,  Northridge,  CA  earthquake  of  January  17,  1994  and  the  S16E  component  of  Pacoima  Dam,  San 
Fernando,  CA  earthquake  of  March  22,  1971  scaled  to  an  EPA  = 0.4g.  Parameters  for  Bouc-Wen 
Model  are:  A = !,/}=/=  0.5,  n = 20,  and  a = 0.02 
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Figure  4.6.  Relative  input  energy  Eir  for  ductility  ratios  //=  1,  2,  3,  and  5,  and  a damping  ratio  £ = 
5%  for  three  accelerograms  with  short  (a  - c),  and  three  accelerograms  with  long  duration  of  strong 
motion  (d  - 0 in  Table  4.1.  Parameters  for  Bouc-Wen  Model  are:  A = 1,  /?=  y=  0.5,  n = 20,  and  a = 
0.02 
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energy  for  Pacoima  Darn  and  Mexico  City  accelerograms  (Figures  4.6.c  and  4.6.f),  it 
generally  overestimates  the  inelastic  input  energy. 

Figure  4.7  shows  that  the  relative  input  energy  for  a ductility  ratio  fi  = 3,  a post-  to  pre- 
yield stiffness  ratio  a = 0.02,  and  damping  ratios  £ = 0,  2,  5,  10,  20,  and  40%  for  three 
accelerograms  with  short  duration  of  strong  motion  (a  - c),  and  three  with  long  duration 
(d  - f)  in  Table  4.1.  The  reason  for  including  damping  ratios  larger  than  5%  is  that  large 
damping  ratios  may  are  encountered  when  using  supplemental  damping  devices  in  a 
structure.  Figure  4.7  shows  that  for  damping  ratios  smaller  than  5%,  damping  has  little 
effect  on  the  input  energy.  For  damping  ratios  greater  than  5%,  on  the  other  hand, 
damping  has  a significant  influence  on  the  input  energy  spectra,  particularly  for  periods 
close  to  Tpe.  A damping  ratio  % = 40%  may  reduce  the  input  energy  by  approximately 
50%  for  periods  close  to  Tpe-  For  very  long  periods  as  the  damping  increases,  the  input 
energy  increases. 

Figure  4.8  shows  the  energy  ratios  for  a ductility  ratio  n = 3 and  damping  ratios  0,  2,  5, 
10,  20,  and  40%  subjected  to  the  360°  component  of  the  1994  Sylmar  record  with  t = 5 
s.  It  is  observed  that  as  the  damping  increases,  the  EhJEnm  and  ( EblEir)m  ratios  decrease, 
Figures  4.8.c  and  4.8.d.  Figure  4.8.c  indicates  that  for  a structure  with  a 5%  damping  ratio 
approximately  75%  of  the  input  energy  is  dissipated  through  hysteretic  action,  and  for  a 
structure  with  a 40%  damping  only  20%  of  the  input  energy  is  dissipated  through 
hysteretic  action.  The  figure  indicates  that  damping  significantly  influences  the 
distribution  of  input  energy  among  the  energy  components. 

Fajfar  et  al.  (1989),  and  Uang  and  Bertero  (1990)  recommended  that  the  input  energy  be 
used  as  a damage  index.  Results  of  the  study  by  Bruneau  and  Wang  (1996)  and  this 
study  (see  Figure  4.7)  indicate  that  damping  ratios  smaller  than  5%  have  a minor 
influence  on  the  input  energy.  On  the  other  hand,  damping  affects  the  damage  potential 
of  structures  to  a large  degree  (see  Figures  4.8.c  and  4.8.d).  Therefore,  the  input  energy 
should  not  be  used  as  a sole  damage  index.  Earthquake  input  energy,  however,  can  be 
used  as  a measure  of  the  intensity  of  ground  motion.  Peak  ground  acceleration  PGA, 
effective  peak  acceleration  EPA  defined  by  Equation  (1.1),  and  Spectrum  Intensity  SI 
defined  by  equation  (1.2)  have  all  been  used  in  measuring  the  intensity  of  ground  motion, 
but  they  do  not  account  for  the  influence  of  the  duration.  Earthquake  input  energy,  on  the 
other  hand,  accounts  for  the  duration  of  strong  motion,  frequency  content,  and  structural 
properties  and  is  more  meaningful  than  PGA , EPA,  and  SI  for  measuring  the  intensity. 
The  maximum  ordinate  of  the  elastic  input  energy  spectrum  with  5%  damping  may  be 
used  for  scaling  the  records  to  the  same  intensity. 

The  hysteretic  behavior  of  a structure  depends  on  the  material  and  the  structural  system 
such  as  moment  resisting  frame,  braced  frame,  etc.  Two  cases  of  material  non-linearities 
are  considered  in  this  study:  (a)  bilinear  behavior  with  A = 1,  J3  = y~  0.5,  and  n = 20,  and 
a = 0.02  proposed  by  Wen  (1976)  for  a moment  resisting  steel  frame  (MRF),  and  (b) 
hysteretic  behavior  with  A = 1.4704,  ft  = 0.042,  y = -0.121,  n = 2,  and  a = 0.2325 
proposed  by  Loh  and  Chung  (1993)  for  an  eccentrically  braced  steel  frame  (EBF).  Figure 
4.9  shows  the  hysteresis  loops  for  an  MRF  and  an  EBF  structure  with  T„  = 2.0  s,  ju  = 3, 
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Figure  4.7.  Relative  input  energy  Eir  for  a ductility  ratio  of  3,  and  damping  ratios  %=  0,  2,  5,  10,  20, 
and  40%  for  three  accelerograms  with  short  duration  of  strong  motion  (a  - c),  and  three 
accelerograms  with  long  duration  (d  - f)  in  Table  4.1.  Parameters  for  Bouc-Wen  Model  are:  .4  = 1,/? 
= y=  0.5,  n = 20,  and  a - 0.02 
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Figure  4.8.  Energy  ratio  spectra  for  bilinear  behavior  with  a ductility  ratio  3 and  damping  ratios 
0%,  2%,  5%,  10%,  20%,  and  40%  subjected  to  the  360°  component  of  Sylmar  County  Hospital 
Parking  Lot  Station,  Northridge  earthquake  of  January  17,  1994.  (a)  Ekrm/Eirm,  (b)  Esm/Eirm,  (c) 
Ehn/Eirm,  (d)  (E,/Eir)m,  (e)  EdJEirm , and  (f)  Edism/Eirm.  Parameters  for  Bouc-Wen  Model  are:  A = = 

y=  0.5,  it  = 20,  and  a = 0.02 
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Figure  4.9.  Force-displacement  diagram  for  an  SDOF  MRF  and  an  SDOF  EBF  with  T„  = 2.0  s,  a 
ductility  ratio  of  3,  and  a 5%  damping  ratio  subjected  to  the  360°  component  of  Sylmar  County 
Hospital  Parking  Lot,  Northridge,  CA  earthquake  of  January  17,  1994.  Parameters  for  Bouc-Wen 
Model  are:  (a )A  = 1 , /?  = y=  0.5,  n = 20,  and  a = 0.02,  and  (b)  A = 1.4704,/?=  0.042,  y=  -0.121,  /;  =2, 
and  a=  0.2325 
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and  a 5%  damping  ratio  subjected  to  the  360°  component  of  the  Sylmar  record.  The  total 
input  energy  transmitted  by  that  record  to  the  MRF  and  EBF  structures  at  a period  of  2 s 
are  approximately  1.9  and  4.5  m2/s2  (see  Figure  4.10) 

Figure  4.10  shows  the  input  energy  spectra  for  a ductility  ratio  of  5,  a 5%  damping  ratio, 
and  the  two  types  of  hysteretic  behavior  for  the  360  component  of  the  Sylmar  record  and 
the  S00E  component  of  the  Mexico  City  record.  The  figure  indicates  that  the  shape  of 
the  hysteresis  loop  has  a significant  influence  on  the  magnitude  and  variation  of  the  input 
energy  spectra. 

The  results  of  this  study  indicate  that  the  influences  of  structural  properties  on  earthquake 
input  energy  are  as  important  as  the  influences  of  ground  motion  characteristics. 
Nonetheless,  the  empirical  equations  proposed  by  Akiyama  (1985),  Kuwamura  and 
Galambos  (1989),  Fajfar  et  al.  (1989),  Uang  and  Bertero  (1990),  and  Kuwamura  et  al. 
(1993)  for  estimating  the  earthquake  input  energy  are  based  on  the  ground  motion 
characteristics  only.  The  equation  proposed  by  Housner  (1956)  does  not  account  for  the 
influence  of  the  duration  of  strong  motion.  The  procedure  proposed  by  Manfredi  (2001) 
for  the  constant-velocity  region  of  response  spectra  accounts  for  the  influences  of  the 
ground  motion  characteristics  as  well  as  structural  properties  in  estimating  the  earthquake 
input  energy. 

A procedure  for  computing  the  earthquake  input  energy  from  the  response  spectra  by 
accounting  for  the  soil  type  is  needed.  A large  number  of  accelerograms  recorded  on 
different  soil  should  be  used  to  compute  the  energy  spectra  for  various  hysteretic 
behaviors. 
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Figure  4.10.  Relative  input  energy  Eir  for  a ductility  ratio  of  3,  damping  ratio  £ = 5%,  and  two 
hysteretic  behaviors  for  (a)  the  360°  component  of  Sylmar  County  Hospital  Parking  Lot,  Northridge, 
CA  earthquake  of  January  17,  1994,  and  (b)  the  S00E  component  of  Mexico  City  earthquake  of 
September  19,  1985.  Parameters  for  Bouc-Wen  Model  are  (I)  for  MRF  structure:  A = 1,  /}=  y=  0.5, 
n = 20,  and  a-  0.02,  and  (II)  for  EBF  structure  A = 1.4704,  f}=  0.042,  y=  -0.121,  n =2,  and  a=  0.2325 
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5. 


HYSTERETIC  ENERGY  DEMAND 


In  a seismic  design  procedure  based  on  energy,  one  should  limit  the  structural  damage  by 
providing  sufficient  ductility  and  energy  dissipation  capacity  by  hysteretic  action  and/or 
damping  in  the  structure.  The  damage  potential  is  associated  with  the  maximum 
hysteretic  energy  demand  during  the  excitation  and  during  the  largest  yield  excursion  (see 
Figures  3.1  - 3.3).  The  hysteretic  energy  demand  can  be  computed  from  the  input  energy 
spectra  if  the  ratio  of  the  maximum  hysteretic  energy  to  the  maximum  input  energy 
EhJEjm,  is  known.  To  examine  the  relationship  between  hysteretic  energy  and  damage 
potential  three  energy  ratios  are  considered:  1)  the  maximum  ratio  of  hysteretic  to  input 
energy  (Ef/Eir)m  generally  occurring  during  the  largest  yield  excursion,  2)  the  ratio  of  the 
maximum  hysteretic  energy  to  the  maximum  input  energy  EhJEirm  occurring  at  the  end 
of  the  excitation,  and  3)  the  equivalent  number  of  yield  excursions  Neq , see  equation 
(4.13). 

Figure  5.1  shows  ( Eh/Eir)m , EhJEirm , and  Neq  for  a ductility  ratio  /j  = 3 and  a 5%  damping 
ratio  for  the  records  with  short  and  long  duration  of  strong  motion.  Figures  5.1. a to  5.1.d 
indicate  that,  on  the  average,  (Eh/Eir)m  are  approximately  20%  larger  than  EhJEirm 
regardless  of  the  period  of  structure.  Ground  motion  characteristics  have  a minor 
influence  on  the  ( Eh!Eir)m  ratio  (Figures  5.1. a and  5.1.b)  and  a major  influence  on  the 
(Eh/Ejr)m  ratio  (Figures  5.1.c  and  5.1.d).  For  example,  for  an  SDOF  structure  with  T„  = 
1.0  s (Figure  5.2)  subjected  to  the  360°  component  of  1994  Sylmar  record  (tj,  = 5 s), 
approximately  90%  of  the  total  hysteretic  energy  is  dissipated  during  the  first  (largest) 
yield  excursion  where  ( Eh/Eir)m  = 0.77  and  EhJEirm  = 0.67.  Similarly  for  a structure  with 
T„  = 2.0  s (Figure  5.3),  70%  of  Ef,m  is  dissipated  in  the  first  (largest)  yield  excursion 
where  (Eh/Etr)m  = 0.71  and  EhmIEirm  = 0.69.  Therefore,  when  a flexible  structure  (long 
period)  is  subjected  to  a short  duration  earthquake,  the  damage  potential  during  the 
largest  yield  excursion,  indicated  by  (Ef,/Eir)m,  may  be  more  important  than  the  damage 
potential  due  to  the  total  number  of  yield  excursions  indicated  by  EhJEirnl. 

Figures  5.1.e  and  5.1.f  show  that  generally  for  the  same  ductility  and  damping,  the 
equivalent  number  of  yield  excursions  Neq  is  larger  for  the  records  with  long  duration  of 
strong  motion  than  for  the  records  with  short  duration,  indicating  a larger  damage 
potential  for  the  long  duration  records.  For  the  two  longest  duration  records,  the  1985 
Mexico  City  - S00E  (/j,  = 71  s)  and  1985  Llolleo  - 100°  (fa  = 41  s),  Neq  is  the  largest  in 
the  vicinity  of  the  predominant  period  where  Tpe  = 2.0  s and  0.52  s for  the  Mexico  City 
and  Llolleo  records,  respectively  (see  Figure  5.1.f).  Figure  5.1  indicates  that  the 
structures  with  T„  =1.0  to  2.5  s and  T„  = 0.2  to  0.7  s have  the  potential  for  significant 
damage  when  subjected  to  the  Mexico  City  and  Llolleo  records,  respectively.  One 
cannot  make  a definitive  conclusion  from  Figures  5.1.e  and  5.1.f,  however,  for  the 
damage  potential  of  structures  subjected  to  the  other  18  accelerograms. 
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Figure  5.1.  (a,  b)  Maximum  ratio  of  the  hysteretic  to  input  energy  (EhIEir)m,  (c,  d)  Ratio  of  the 
maximum  hysteretic  energy  to  the  maximum  input  energy  EhJEirm,  (e,  f)  Equivalent  number  of  yield 
excursions  Neq  for  records  with  short  and  long  durations  for  a ductility  ratio  fi  = 3 and  a 5%  damping 
ratio.  Parameters  for  Bouc-Wen  Model  are:  A = 1,  f}=  /=  0.5,  n = 20,  and  a - 0.02. 
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Figure  5.2.  Time-histories  of  EhIEhm  and  EhIEir  for  a bilinear  SDOF  structure  with  T„  = 1.0  s,  a 
ductility  ratio  //  = 3,  and  5%  damping  subjected  to  the  90°  component  of  Sylmar  Count>  Hospital 
Parking  Lot,  Northridge,  CA  earthquake  of  January  17,  1994.  Parameters  for  Bouc-Wen  Model  are 
A = 1,  /?  = y=  0.5,  n = 20.  a=  0.02 
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Figure  5.3  demonstrates  that,  in  general,  the  energy  ratios  considered  provide  different 
information  about  the  damage  potential  of  a structure.  The  ratio  ( Eh/E,r)m  reflects  the 
damage  due  to  the  largest  yield  excursion  and  is  independent  of  the  duration  of  strong 
motion  and  period  of  structure.  The  energy  ratios  EhJEirm  and  Neq,  on  the  other  hand, 
reflect  the  damage  potential  due  to  cumulative  inelastic  deformation  for  the  entire 
duration  of  accelerogram.  The  equivalent  number  of  yield  excursions  Neq  depends  on  the 
duration  of  strong  motion  and  the  period  of  structure,  whereas  the  EhJEirm  ratio  is 
independent  of  both  only  for  periods  less  than  1 5.  For  quantifying  the  damage  potential, 
the  Ehm/Eim,  and  {E\JEir)m  ratios  may  be  used  to  define  a damage  index  to  address  the 
damage  potential  for  the  largest  yield  excursion  and  the  total  number  of  yield  excursions 
during  the  earthquake. 

Figures  5.4. a and  5.4.b  show  the  mean  (Ey,IEir)m  and  E\mJEwm  ratios,  respectively,  for  the 
20  accelerograms  for  5%  damping  and  ductility  ratios  /u  = 2,  3,  and  5.  In  general,  as  the 
ductility  increases  both  ratios  increase,  indicating  a larger  dissipation  of  input  energy 
through  hysteretic  action.  Figure  5.4.b  can  be  used  to  estimate  the  hysteretic  energy 
spectra  from  input  energy  spectra  if  the  relationship  among  the  Eh,JEirm  ratio,  the  duration 
of  strong  motion,  and  the  period  of  structure  is  determined.  A large  number  of 
accelerograms  recorded  on  different  soil  should  be  used  and  energy  spectra  for  various 
hysteretic  behaviors  should  be  computed  to  develop  a procedure  for  estimating  the 
hysteretic  energy  demand. 
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Figure  5.3.  Time-histories  of  EhIEhm  and  EhIEir  for  a bilinear  SDOF  structure  with  T„  = 2.0  s,  a 
ductility  ratio  n = 3,  and  5%  damping  subjected  to  the  90°  component  of  Sylmar  County  Hospital 
Parking  Lot,  Northridge,  CA  earthquake  of  January  17,  1994.  Parameters  for  Boue-Wen  Model  are 
A = \^p=y-  0.5,  n = 20,  a = 0.02 
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Figure  5.4.  The  mean  ( EhIEir)m  and  EhJEirm  spectra  for  bilinear  behavior  with  5%  damping  and 
ductility  ratios  p = 2,  3,  and  5 for  20  accelerograms.  Parameters  for  Bouc-Wen  Model  are:  A = 1,  /?  = 
y=  0.5,  n = 20,  and  a=  0.02. 
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6. 


ENERGY  RATIOS  IN  AN  MDOF  STRUCTURE 


Studies  by  McKevitt  et  al.  (1980)  and  Nakashima  et  al.  (1996)  suggest  that  the  energy 
transmitted  to  an  MDOF  structure  can  be  estimated  from  the  energy  transmitted  to  an 
equivalent  SDOF  structure.  Once  the  input  and  the  hysteretic  energies  for  a structure  are 
estimated,  one  needs  to  distribute  the  energies  through  the  height  of  the  structure. 
Sufficient  strength  and  energy  dissipation  capacity  should  be  provided  in  structural 
members  to  limit  the  damage  potential  of  the  structure  to  acceptable  thresholds  (refer  to 
Vision  2000).  This  Chapter  examines  the  distribution  of  earthquake  input  energy  and  its 
components  in  a five-story  building  using  the  20  accelerograms  in  Table  4.1  scaled  to  an 
effective  peak  acceleration  EPA  = 0.4o\ 

The  structural  properties  for  the  building  are  as  follows:  all  floors  have  a weight  W = 
1000  KN.  The  pre-yield  stiffness  of  the  first  story  is  Kj  - 228.898  KN/m  and  the 
stiffnesses  for  stories  2-5  are  reduced  linearly  by  10%,  i.e.  K/K]  = 0.9,  0.8,  0.7,  and  0.6 
for  j = 2 to  5,  respectively.  Damping  is  assumed  proportional  to  stiffness.  Considering  a 
5%  damping  in  the  first  mode  of  vibration  of  the  elastic  structure,  the  damping  coefficient 
for  the  first  story  is  C\  - 3643  KN-s/m  and  the  damping  distribution  through  the  height  is 
C/Ci  = 0.9,  0.8,  0.7,  and  0.6  for  j = 2 to  5,  respectively.  The  yield  deformation  for  each 
story  varies  proportionally  to  the  stiffness  (Yang  et  al.  1992),  where  Uyj/uyi  = 0.9,  0.8,  0.7, 
and  0.6  for  j = 2 to  5,  respectively,  and  uyi  is  the  yield  deformation  for  the  first  story.  The 
pre-yield  periods  of  the  building  are  T„  = 0.50,  0.17,  0.11,  0.08,  and  0.07  s.  The  Bouc- 
Wen  parameters  (Appendix  A)  used  in  the  analyses  are  A = 1,  J3=  y=  0.5,  n = 6 , a = 0.02 
for  all  floors.  The  energy  balance  equation  for  MDOF  structures  is  presented  in 
Appendix  B. 

Figure  6.1  shows  the  total  hysteretic  and  input  energies  for  the  five-story  building  with 
uyi  - 0.5  cm  vs.  the  intensity-based  duration  computed  for  the  20  records.  As  expected, 
the  figure  shows  that  the  total  input  and  hysteretic  energies  for  the  records  with  long 
duration  are  greater  than  those  for  the  records  with  short  duration.  The  standard 
deviations  of  Eh  for  the  records  with  short  and  long  duration  are  0.1  MN.m  and  0.4 
MN.m,  respectively,  and  the  standard  deviations  of  E,r  for  the  short  and  long  duration 
records  are  0.1  MN.m  and  0.5  MN.m,  respectively.  The  large  dispersion  from  the  means 
in  Figure  6.1  is  caused  by  the  influence  of  the  frequency  content  on  input  energy,  i.e.  the 
input  and  hysteretic  energies  for  a structure  with  a period  approximately  equal  to  the 
predominant  period  Tpe,  see  Table  4.1,  are  larger  than  those  for  a structure  with  period 
different  than  the  predominant  period.  Figure  6.1  indicates  that  the  duration  of  strong 
motion  has  a significant  influence  on  energy  demand. 

Table  6.1  shows  the  inter-story  ductility,  normalized  story  shear,  energy  distribution,  and 
the  ratios  (Eh/Eir)m,  Ehm/E,rm , and  Edm/Ejrm  for  each  story  when  the  building  is  subjected  to 
the  S00E  component  of  El  Centro  accelerogram  scaled  to  an  EPA  = 0.4g.  The  yield 
deformation  for  the  first  story  is  ux\  - 0.5  cm.  Table  6.1  shows  that  the  ductility  demand 
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Figure  6.1.  Total  hysteretic  and  input  energies  for  the  five-story  building  with  uvI  = 0.5  cm  subjected 
to  the  20  accelerograms  scaled  to  an  EPA  = OAg.  Parameters  for  Bouc-Wen  Model  are  .4  = 1,  J3=  y= 
0.5,  n - 6,  a - 0.02 
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Table  6.1.  Inter-story  ductility,  normalized  story  shear,  percentage  of  energy  in  each  story,  and  the 
energy  ratios  (Ei/Eir)m,  EhmIEirm , and  EdmIEdm  for  the  five-story  building  subjected  to  the  SOOE 
component  of  El  Centro,  the  Imperial  Valley,  CA  earthquake  of  May  18,  1940  scaled  to  an  EPA  = 
0.4 g. 
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is  the  largest  in  the  third  story  (column  2),  whereas  the  hysteretic  energy  demand  is  the 
largest  in  the  first  story  (column  7).  The  ductility  ratio  does  not  consider  the  influence  of 
the  duration,  whereas  the  hysteretic  energy  does.  Therefore,  it  is  reasonable  that  in 
design  both  ductility  and  energy  demands  be  examined  for  each  story.  It  is  observed  that 
for  higher  stories  the  relative  kinetic  energy  Ef  is  larger  (column  4)  since  it  is  computed 
from  the  relative  velocity,  whereas  elastic  strain  energy  Es  is  smaller  (column  6)  since  it 
is  computed  from  the  drift*5  which  for  a shear  building  is  smaller  in  higher  stories.  The 
maxima  of  Ef/Eir  are  presented  in  column  10  of  Table  6.1  which  are  approximately  equal 
for  the  four  lower  stories.  The  fifth  story  behaves  approximately  elasticly  (ju  = 1.3  in 
column  2).  Only  9%  of  the  input  energy  is  dissipated  through  hysteresis  (column  1 1 ) and 
51%  of  the  energy  is  dissipated  through  damping.  It  is  also  observed  that  (Ef/EJ„,  is 
generally  larger  than  EhJE,rm  for  each  story  and  the  difference  between  the  two  ratios 
increases  for  higher  stories  which  have  smaller  ductility  demands. 

Table  6.2  shows  the  influence  of  the  yield  deformation  on  the  maximum  inter-story 
ductility,  base  shear,  and  total  energy  parameters  for  the  five-story  building  subjected  to 
the  S00E  component  of  the  1 940  El  Centro  accelerogram  scaled  to  an  EPA  = OAg.  It  is 
observed  that  as  the  yield  deformation  decreases  ductility  demand  increases  (column  2), 
and  base  shear  decreases  significantly  (column  3),  whereas  the  total  dissipated  and  input 
energies  remain  approximately  the  same  (columns  8 and  9).  Nonetheless,  ductility  has  a 
pronounced  influence  on  the  distribution  of  input  energy  among  energy  components, 
specially  for  ductility  demands  less  than  3 (columns  4 - 7). 

Figure  6.2  shows  the  relative  displacement,  ductility  demand,  normalized  story  shear,  the 
energy  ratio  EhmIEirm,  and  the  distribution  of  the  hysteretic  and  input  energies  though  the 
height  of  the  5-story  building  with  yield  deformations  uvi  = 0.3,  0.5,  0.7,  1.0,  and  1.5  cm, 
subjected  to  the  S00E  component  of  the  1940  El  Centro  accelerogram  scaled  to  an  EPA  = 
OAg.  Also  presented  are  the  results  for  the  elastic  case.  The  figure  indicates  that  as  the 
yield  deformation  decreases  the  relative  displacement  for  each  story  decreases,  and  the 
ductility  and  hysteretic  energy  demands  increase.  It  is  observed  that  structural  non- 
linearity has  a small  effect  on  the  distribution  of  story  shear,  hysteretic  energy,  and  input 
energy  for  this  building,  see  Figures  6.2.c,  6.2.e,  and  6.2.f.  Even  though  displacement  is 
largest  at  the  top,  the  ductility  and  hysteretic  energy  demands  are  largest  in  the  lower 
three  stories. 

Figure  6.3  compares  the  percentage  of  energy  distributed  through  the  height  of  the  five- 
story  building,  with  a yield  deformation  for  the  first  story  uyj  — 0.5  cm,  subjected  to  the 
20  accelerograms  scaled  to  an  EPA  = OAg.  The  figure  indicates  that  in  general,  the 
distribution  of  energy  through  the  height  of  the  building  is  not  affected  by  ground  motion 
characteristics.  On  the  average,  23,  28,  27,  15,  and  7%  of  the  total  input  energy  is 
transmitted  to  stories  1 to  5,  respectively.  Similarly,  22,  30,  31,  16,  and  1%  of  the  total 
hysteretic  energy  is  dissipated  in  stories  1 to  5,  respectively.  While  the  above 
percentages  provide  rough  estimates  and  trend  in  distribution  and  dissipation  of  energy,  a 
large  number  of  structures  and  records  should  be  examined  to  arrive  at  a rational  formula 
for  distribution  of  energy  through  the  height  of  a structure. 

b The  difference  between  two  consecutive  floor  displacements. 
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Table  6.2.  Influence  of  yield  deformation  on  maximum  inter-story  ductility  base  shear  Qb , and 
total  energy  for  the  five-story  building  subjected  to  the  SOOE  component  of  El  Centro,  the  Imperial 
Valley,  CA  earthquake  of  May  18, 1940  scaled  to  an  EPA  = 0.4g. 
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Figure  6.2.  (a)  Relative  displacement  ureh  (b)  ductility  demand  //,  (c)  normalized  story  shear,  (d) 
Ehn/Eirm , (e)  distribution  of  Eh  though  the  height,  and  (0  distribution  of  Eir  through  the  height  for  the 
5-story  building  with  yield  deformations  of  the  first  story  uyl  = 0.3,  0.5,  0.7,  1.0,  and  1.5  cm  subjected 
to  S00E  component  of  El  Centro,  the  Imperial  Valley,  CA  earthquake  of  May  18,  1940  scaled  to  an 
EPA  = 0.4g.  Also  presented  are  the  results  for  elastic  behavior.  Parameters  for  Bouc-Wen  Model 
are:  A = 1,  f}=  y=  0.5,  n = 20,  and  a - 0.02 
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Figure  6.3  Distribution  of  energy  through  the  height  of  a 5-story  building  with  uyl  = 0.5  cm  subjected 
to  20  accelerograms  scaled  to  an  EPA  = 0.4 g.  Parameters  for  Bouc-Wen  Model  are:  A = 1,  f}=  y = 
0.5,  /;  = 20,  and  a=  0.02 
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7.  DISTRIBUTION  OF  INPUT  ENERGY  IN  STRUCTURES  WITH 
BASE-ISOLATION,  SUPPLEMENTAL  DAMPING  AND  SEMI- 
ACTIVE CONTROL 


Base  isolation.  Figure  7.1.b,  supplemental  damping,  and  semi-active  control  are  used  to 
reduce  the  damaging  effects  of  earthquakes.  Base  isolation  reduces  the  motion 
transmitted  to  the  structure  and  the  base  shear  by  increasing  the  period  of  structure. 
Supplemental  damping  reduces  the  hysteretic  energy  demand  on  the  structure  by 
dissipating  a portion  of  the  input  energy.  Semi-active  control  combines  the  features  of 
active  and  passive  control.  Similar  to  passive  systems,  control  forces  are  generated  by 
the  response  of  the  building,  and  similar  to  active  systems,  controllers  are  used  to  monitor 
feedbacks  and  apply  the  control  forces.  Unlike  active  systems  which  require  a large 
power  source  to  generate  the  control  forces,  semi-active  systems  require  a smaller  power 
source. 

Semi-active  control  (Figure  7.1.c)  includes  active  variable  stiffness  and  active  variable 
damping.  In  variable  stiffness  control,  the  stiffness  of  the  building  is  adjusted  to  prevent 
resonance  of  the  building  from  the  excitation.  In  variable  damping  control,  supplemental 
energy  dissipation  devices  such  as  friction  and  fluid  dampers  are  used  to  reduce  the 
structural  response.  This  chapter  presents  the  energy  balance  equation  for  structures  with 
base-isolation  and  semi-active  control.  Section  7.1,  and  the  energy  distribution  for  four 
five-story  buildings  with:  fixed-base,  base-isolation,  supplemental  damping,  and  semi- 
active control  in  Section  7.2. 


7.1.  Energy  Equations  for  Buildings  with  Base-Isolation  and  Semi- 
Active  Control 

The  energy  balance  equation  for  a base-isolated  building  with  a base  mass  nih  (Figure 

7.1. b)  is  similar  to  a fixed-base  building  (Figure  7.1. a)  with  an  additional  degree  of 
freedom  for  the  base  isolation,  see  equation  B.2  in  Appendix  B.  Assuming  that  the 
building  undergoes  a rigid  body  motion,  the  period  of  the  base-isolated  building  can  be 
approximated  as 


+ mb 

and 


(7.1) 


and  the  damping  ratio  of  the  base  isolation  system  as 
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Figure  7.1.  (a)  Five-story  building  with  fixed  base,  (b)  with  base-isolation,  and  (c)  with  supplemental 
damping  or  semi-active  control  subjected  to  the  S00E  component  of  El  Centro,  the  Imperial  Valley, 
CA  earthquake  of  May  18,  1940  scaled  to  an  EPA  = 0.4g 
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(7.2) 


5 


The  energy  balance  equation  for  a building  with  supplemental  damping  (Figure  7.1.c)  is 
similar  to  a fixed-base  building,  see  equation  B.2  in  Appendix  B.  The  total  damping  of 
the  structure  (inherent  and  supplemental)  should  be  used  in  the  analysis. 

The  governing  differential  equation  of  motion  for  a non-linear  »-story  shear  building  with 
m semi-active  dampers  (Figure  7.1.c)  is  (see  Yang  et  al.  1992) 


Ft  (0  + Fd  (t)  + Fs  (0  = ug  ( t ) + DFC  (t) 


(7.3) 


where  Ft  is  an  n - dimensional  inertia  force  vector,  Fd  is  an  Yl  - dimensional  damping 
force  vector,  Fs  is  an  Yl  - dimensional  non-linear  restoring  force  vector,  is  an  Yl  - 


dimensional  mass  vector,  ii  is  ground  acceleration,  D is  an  Yixm  matrix  defining  the 


locations  of  control  forces,  and  Fc  denotes  an  m -dimensional  control  force  vector.  The 
inertia  force  vector  can  be  written  as  Ff  = M ii  where  M is  an  nx  n mass  matrix  and 


ii  is  an  Yl-  dimensional  acceleration  vector.  The  damping  force  can  be  written  as 


Fd  — C it  where  C is  an  n x n damping  matrix  and  ii  is  an  Yl  - dimensional  velocity 
vector.  The  restoring  force,  however,  is  a non-linear  function  of  u and  ii . 


Using  the  state-space  representation,  equation  (7.3)  can  be  written  as 


Of)  = g[z(?)]  + B u(t)  + H ug  (0 


(7.4) 


in  which  the  2yi  - dimensional  vector  g[z(t )]  is  given  by 


(7.5) 


and  z{t ) is 


(7.6) 


The  control  input  matrix  B is 
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(7.7) 


0 

M lD 


and  the  disturbance  input  vector  H is 


H = 


(7.8) 


For  the  linear  quadratic  regulator  algorithm  (LQR),  the  optimal  control  force  Fc(t ) in 
equation  (7.3)  is  selected  by  minimizing  the  performance  index  J 

J=^\zT(t)Qz(t)  + Fj(t)RFc(t)]dt  (7.9) 


where  t f is  the  duration  of  accelerogram,  Q — ql  is  a 2/7  x 2/7  symmetric  positive 

semi-definite  weighting  matrix  ( q is  a scalar  and  / is  unit  matrix),  and  R = rl  is  an 

777  x 777  positive  definite  weighting  matrix  ( r is  a scalar  and  / is  unit  matrix).  By 
linearizing  the  hysteretic  behavior  of  the  structure  at  the  initial  equilibrium  state  and 
minimizing  the  performance  index  subject  to  the  constraint  of  equation  (7.4),  the  control 
forces  are  obtained  as 

Fc (0  = -0.5 R~'BrPz{t)  = G :(t)  (7.10) 


where  G is  an  777  x 2/7  feedback  gain  matrix  and  P is  a 2/7  x 2/7  matrix  obtained  from 
the  solution  of  the  classical  Riccati  equation 

2Q  z + A(z)P  z + P z + P z = 0 (7.11) 

where 

A (z)  = dg(z)ldz  (7.12) 

Transposing  equation  (7.3),  multiplying  both  sides  by  dll  = il  dt  and  integrating,  one 
obtains  the  energy  balance  equation  as 

JV,  (t)du+  jUQF/(t)du+  jy/(t)du  = -jjug(t)du+jUQFcT(<t)DT  da  (7.13) 

The  energy  corresponding  to  the  controller  is 
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Ec  = JV/  DTi,  dt 


(7.14) 


Other  energy  vectors  are  defined  in  Appendix  B.  Equation  (7.13)  can  be  written  as 

Ekr  + Ed  + Es  + Eh  ~ Eir  + Ec  (7-15) 

where  each  term  is  an  n -dimensional  vector. 


7.2.  Numerical  results 

The  lead-core  rubber  isolation  bearing  used  in  the  example  has  a period  Tb  = 2.0  s and  a 
damping  ratio  %b  = 10%.  The  weight  of  the  base  is  equal  to  the  weight  of  individual 
floors  ft^=1000  KN,  the  stiffness  of  the  base  isolation  system  kb  = 6036  KN/m,  and  the 

damping  coefficient  cb  = 384  KN-s/m.  The  Bouc-Wen  parameters  for  the  lead-core 

rubber  bearing  are  A = 1,  /?  = y=  0.5,  n = 3,  and  a = 0.6  (Yang  et  al.,  1992).  For  the 
building  with  supplemental  damping  one  damper  is  used  at  each  story  and  the  total 
damping  (inherent  and  supplemental)  is  40%.  For  the  building  with  semi-active  control 
one  damper  is  used  at  each  story.  The  weight  parameters  used  for  20  accelerograms  are  r 
= 10'5  and  q = 10+5.  As  q increases,  the  control  force  increases.  The  yield  deformation 
for  the  first  story  is  uv]  = 0.5  cm  for  all  buildings. 

Figure  7.2,  7.3,  and  7.4  show  the  percentage  of  energy  distributed  through  the  height  of 
the  five-story  building  with  base  isolation,  supplemental  damping,  and  semi-active 
control  subjected  to  the  20  records  scaled  to  an  EPA  = 0.4g.  Figures  6.3,  7.2,  7.3,  and  7.4 
conclude  that  ground  motion  characteristics  have  a minor  influence  on  the  distribution  of 
energy  through  the  height  of  the  buildings.  Figure  7.2  shows  that  base  isolation  is 
effective  for  structures  subjected  to  all  accelerograms,  except  the  Mexico  City 
accelerogram  since  the  period  of  base-isolation  {Tb  = 2.0  s)  is  close  to  the  predominant 
period  of  the  accelerogram  (Tpe  - 2.0  s).  Therefore,  using  base-isolation  for  a 5 story 
building  in  Mexico  City  is  detrimental,  while  it  is  effective  in  reducing  the  ductility  and 
energy  dissipation  demand  for  building  in  other  locations. 

The  average  percentage  of  hysteretic  energy,  input  energy,  and  story  shear  in  each  story 
for  the  four  buildings  subjected  to  the  20  accelerograms  scaled  to  an  EPA  - OAg  are 
given  in  Table  7.1  and  plotted  in  Figure  7.5.  The  figure  indicates  that  the  energy 
distribution  for  the  buildings  with  fixed  base,  supplemental  damping,  and  semi-active 
control  are  similar,  whereas  the  energy  distribution  for  the  building  with  base  isolation  is 
significantly  affected  by  the  isolator,  i.e.  the  majority  of  the  input  energy  is  absorbed  and 
dissipated  in  the  isolator.  The  story  shear  distribution  is  also  similar  for  the  four 
buildings. 

Table  7. 2. a presents  the  maximum  drift,  the  maximum  relative  displacement,  the  ductility 
ratio,  and  the  maximum  story  shear,  and  Table  7.2.b  shows  the  maximum  energy 
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5 


Figure  7.2  Distribution  of  energy  through  the  height  of  a 5-story  building  with  base  isolation  (10% 
damping,  Tb  = 2.0  s,  and  urb  = 4.0  cm)  with  uyI  = 0.5  cm  subjected  to  20  records  scaled  to  an  EPA  = 
0.4g.  Parameters  for  Bouc-Wen  Model  for  the  stories  are:  A = 1,  /?=  y=  0.5,  n = 20,  and  a=  0.02  and 
for  base  isolation  system  are:  A - 1 , /?  = y = 0.5,  n - 3,  and  a = 0.6 
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Figure  7.3.  Distribution  of  energy  through  the  height  of  a 5-story  building  with  supplemental 
damping  (total  damping  40%)  with  uyl  = 0.5  cm  subjected  to  20  records  scaled  to  an  EPA  = 0.4g. 
Parameters  for  Bouc-Wen  Model  for  the  stories  are:  A = 1,  /?=  y=  0.5,  n = 20,  and  a - 0.02 
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Figure  7.4.  Distribution  of  energy  through  the  height  of  a 5-story  building  with  semi-active  control  (r 
= 10'5  and  q = 10+5)  with  uyI  = 0.5  cm  subjected  to  20  records  scaled  to  an  EPA  = 0.4#.  Parameters  for 
Bouc-Wen  Model  for  the  stories  are:  A = 1,  /?=  y=  0.5,  n = 20,  and  a=  0.02 
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Table  7.1.  The  average  distribution  percentage  of  hysteretic  energy,  input  energy,  and  story  shear,  in 
percent,  in  each  story  of  the  five-story  building  with  fixed-base,  base  isolation,  supplemental 
damping,  and  semi-active  control  subjected  to  20  accelerograms  scaled  to  an  EPA  = 0.4g. 


Tabl^J^Distributioi^ni^stereti^nerg^£^%) 


1 

2 

3 

4 

5 

Story 

Fixed 

Base 

Suppl. 

Semi 

Base 

Isolation  Damp. 

Active 

5 

1 

0 

1 

1 

4 

16 

0 

8 

12 

3 

31 

0 

22 

30 

2 

30 

1 

33 

34 

1 

22 

1 

36 

23 

Base 

- 

98 

- 

- 

Table  7.1. b Distribution  of  in 

put  energy 

Eir(%) 

1 

2 

3 

4 

5 

Story 

Fixed 

Base 

Suppl. 

Semi 

Base 

Isolation 

Damp. 

Active 

5 

7 

6 

6 

7 

4 

15 

6 

11 

20 

3 

27 

6 

20 

21 

2 

28 

6 

28 

19 

1 

23 

6 

35 

23 

Base 

- 

71 

- 

- 

Table  7.1.c  Distribution  of  story  shear 

(?(%) 

1 

2 

3* 

4 

5 

Story 

Fixed 

Base 

Suppl. 

Semi 

Base 

Isolation 

Damp. 

Active 

5 

32 

26 

29 

31 

4 

50 

49 

48 

49 

3 

67 

68 

64 

66 

2 

84 

85 

82 

83 

1 

100 

100 

100 

100 

The  distribution  of  stoiy  shear  for  the  building  with  base  isolation  is  given  for  the  super  structure. 
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Figure  7.5.  Distribution  of  (a)  hysteretic  energy  Eh , (b)  input  energy  Eir , and  (c)  story  shear  Q 
distributed  in  the  stories  of  the  5-story  building  with  fixed  base,  base  isolation,  supplemental 
damping,  and  semi-active  control  subjected  to  20  records  scaled  to  an  EPA  = 0.4g.  Parameters  for 
Bouc-Wen  Model  for  the  stories  are:  A = 1,  fi  = y = 0.5,  n = 20,  and  a - 0.02.  Note  that  for  the 
building  with  base  isolation  the  percentage  is  given  for  the  base  shear  of  the  super  structure. 
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Table  7.2.  (a)  Maximum  drift,  maximum  relative  displacement,  ductility  ratio,  and  maximum  story  shear  for  each  story,  and  (b)  maximum  damping, 
maximum  hysteretic,  maximum  controller,  and  maximum  input  energies  for  the  four  buildings  (with  fixed-base,  base-isolation,  supplemental  damping, 
and  semi-active  control)  subjected  to  the  SOOE  component  of  El  Centro,  the  Imperial  Valley,  CA  earthquake  of  May  18, 1940  scaled  to  an  EPA  = 0.4g. 
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Total  input  energy  transmitted  to  stories  1 to  5 are  103  KN-m. 


parameters:  damping,  hysteretic,  controller,  and  input  for  each  story  of  the  four  buildings 
subjected  to  the  SOOE  component  of  El  Centro,  the  Imperial  Valley,  CA  earthquake  of 
May  18,  1940  scaled  to  an  EPA  - OAg.  To  have  a larger  force  in  controllers,  the  weight 
parameters  used  in  the  analysis  are  selected  as  r = 10'5  and  q = 10+6.  Except  for  the  base- 
isolated  building,  the  total  input  energy  is  equal  to  the  sum  of  the  maximum  input 
energies  of  all  stories.  For  the  base-isolated  building,  all  stories  behave  elastically,  and 
the  base  isolation  system  behaves  inelastically.  The  maximum  input  energy  for  the 
stories  occurs  in  the  initial  stages  of  the  excitation,  whereas  the  maximum  input  energy 
for  the  base  isolation  system  occurs  at  the  end  of  the  excitation.  Consequently,  one 
cannot  add  up  the  input  energies  measured  at  different  times  to  obtain  the  total  input 
energy.  The  total  input  energy  for  stories  1 to  5 of  the  base-isolated  building  is  1 03  KN- 
m (see  Table  7.2. b).  The  total  input  energy  transmitted  to  the  buildings  with  base- 
isolation,  supplemental  damping,  and  semi-active  control  are  approximately  20%,  65%, 
and  25%  of  the  total  input  energy  transmitted  to  the  building  with  fixed  base, 
respectively. 

For  the  base-isolated  building  ( Tb  - 2.0  s)  subjected  to  the  El  Centro  accelerogram  (Tpa  = 
0.26  s and  Tpe  = 0.86  s,  see  Table  4.1),  the  base  shear  (see  columns  5 and  9 of  Table  7. 2. a 
and  Figure  7.6.c),  and  the  total  input  energy  (columns  4 and  7 of  Table  7.2.b)  are 
approximately  60%  of  those  for  the  fixed  base  building.  The  story  drifts  (columns  2 and 
6 of  Table  7. 2. a)  of  the  base-isolated  building  are  less  than  30%  of  the  story  drift  for  the 
fixed-base  building.  The  base-isolated  building  behaves  like  a rigid  body  (Figure  7. 6. a) 
and  all  stories  remain  elastic  during  the  excitation  (column  8 of  Table  7. 2. a;  Figures 
7.6.b,  7.6. d,  and  7.6.e).  The  results  indicate  that  the  input  energy  is  dissipated  in  the  base 
isolation  system  (columns  5 and  7 of  Table  7.2.b). 

The  base  shear  for  the  building  with  semi-active  control  is  approximately  90  percent  of 
the  base  shear  for  the  fixed-base  building,  whereas  the  input  energy  is  approximately  20 
percent.  The  maximum  control  forces  generated  in  the  buildings  with  semi-active  control 
is  288  KN. 

Table  7.3  shows  the  distribution  of  energy  among  various  energy  components  for  the  four 
buildings  subjected  to  the  SOOE  component  of  El  Centro,  the  Imperial  Valley,  CA 
earthquake  of  May  18,  1940  scaled  to  an  EPA  = OAg.  The  results  indicate  that  providing 
base  isolation,  supplemental  damping,  or  semi  active  dampers  has  a significant  influence 
on  the  energy  demand. 

Finally,  the  first  story  drift  and  total  input  energy  for  the  four  buildings  subjected  to  the 
SOOE  component  of  the  1940  El  Centro  accelerogram  scaled  to  an  EPA  = OAg  are 
compared  in  Figures  7.7  and  7.8,  respectively.  It  is  observed  that  both  the  drift  and  input 
energy  of  the  buildings  with  base-isolation,  supplemental  damping,  and  semi-active 
control  are  less  than  those  for  the  fixed-base  building.  According  to  Figure  7.7,  base- 
isolation  is  the  most  effective  system  in  reducing  the  base  shear  and  drift  for  this  building 
subjected  to  El  Centro  accelerogram.  Results  conclude  that  base-isolation,  supplemental 
damping,  and  semi-active  control  are  all  effective  in  reducing  the  hysteretic  energy 
demand  for  this  building  subjected  to  El  Centro  accelerogram. 
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Figure  7.6.  (a)  Relative  displacement,  (b)  ductility  demand,  (c)  story  shear,  (d)  Ehn/E,rm , (e)  ( Ei/Eir)m , 
and  (1)  Edn/Eirm  for  the  five-story  buildings:  fixed-base,  base-isolated,  semi-active  with  large  control 
force,  and  semi-active  with  small  control  force  with  uyI  = 0.5  cm  and  = 5%  subjected  to  the  S00E 
component  of  El  Centro,  the  Imperial  Valley  earthquake  of  May  18,  1940  scaled  to  an  EPA  = 0.4g. 
Parameters  for  Bouc-Wen  Model  are:  A = 1,  fi  = y = 0.5,  n = 20,  and  a - 0.02,  and  for  the  base 
isolation  system  are:  A = 1,  /?=  y=  0.5,  n = 3,  and  a - 0.6 
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Table  7.3.  The  energy  ratios  EkrIEir,  EdIEin  EJEir , Eh/Eir , and  EJEir  for  the  the  four  buildings  (a) 
fixed-base,  (b)  base-isolated,  (c)  semi-active  with  small  control  forces,  (d)  and  semi-active  with  large 
control  forces  subjected  to  the  SOOE  component  of  El  Centro,  the  Imperial  Valley,  CA  earthquake  of 
May  18, 1940  scaled  to  an  EPA  = 0.4g. 


Tabl^j3j^^h^uildin^wiidiJlxed-base 


1 

2 

3 

4 

5 

6 

Story 

EkrJEirni 

EsJEirm 

EtlnJ  Ejrm 

Ehm/Eirni 

Edism/Eirm 

1 

0.90 

0.51 

0.03 

0.09 

0.60 

2 

0.26 

0.43 

0.02 

0.57 

1.00 

3 

0.09 

0.33 

0.01 

0.67 

1.00 

4 

0.04 

0.35 

0.02 

0.65 

1.00 

5 

0.01 

0.42 

0.02 

0.58 

1.00 

Total 

0.10 

0.38 

0.02 

0.62 

1.00 

Table  7.3. b The  building  with  base  isolation  system 

1 

2 

3 

4 

5 

6 

Story 

E krnJ  E jrni 

E sm!  E jrrn 

EjJ  E,rw 

E/,JEirm 

E disnJ E jrw 

1 

1.00 

0.00 

0.00 

0.00 

0.00 

2 

0.99 

0.01 

0.01 

0.00 

0.01 

3 

0.98 

0.02 

0.03 

0.00 

0.02 

4 

0.97 

0.02 

0.04 

0.00 

0.02 

5 

0.96 

0.03 

0.06 

0.00 

0.03 

Base 

0.05 

0.71 

0.15 

0.29 

1.00 

Total 

0.33 

0.71 

0.16 

0.29 

1.00 

Table  7.3.c  The  bui 

ding  with  su 

pplemental  damping 

1 

2 

3 

4 

5 

6 

Story' 

EkrJEirm 

EsJEirm 

E(lnJ  Eirm 

E hnJ  E ir  m 

E disnJ  Eirm 

1 

0.32 

0.93 

0.04 

0.02 

0.95 

2 

0.08 

0.93 

0.02 

0.07 

1.00 

3 

0.03 

0.90 

0.02 

0.10 

1.00 

4 

0.01 

0.89 

0.02 

0.11 

1.00 

5 

0.00 

0.90 

0.02 

0.10 

1.00 

Total 

0.03 

0.90 

0.02 

0.10 

1.00 

Table  7.3.d  The  building  with  semi-active  control 


1 

2 

4 

3 

5 

6 

7 

Story 

E krnJ  E jrni 

EsJEirm 

EdmlEjm, 

Eh  J Eirm 

EdisJEirm 

E cm ! E irm 

1 

0.37 

0.12 

0.03 

0.00 

0.95 

0.82 

2 

0.16 

0.21 

0.05 

0.07 

1.00 

0.72 

3 

0.06 

0.25 

0.05 

0.19 

1.00 

0.56 

4 

0.02 

0.28 

0.06 

0.22 

1.00 

0.50 

5 

0.01 

0.32 

0.08 

0.14 

1.00 

0.54 

Total 

0.08 

0.26 

0.06 

0.15 

1.00 

0.59 
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Figure  7.7.  Drift  time-histories  for  the  first  story  of  the  buildings  with  fixed-base,  base-isolation, 
supplemental  damping,  and  semi-active  control  subjected  to  the  SOOE  component  of  El  Centro,  the 
Imperial  Valley,  CA  earthquake  of  May  18,  1940  scaled  to  an  EPA  = 0.4g.  Parameters  used  for 
Bouc-Wen  Model  for  the  stories  are  A = 1,  f}=  y—  0.5,  n = 20,  a = 0.02,  and  for  the  base  isolation 
system  are  A = 1,  /?=  y=  0.5,  n = 3,  a=  0.6.  The  total  damping  of  the  building  with  supplemental 
damping  is  40%.  The  weighting  matrices  for  semi-active  building  are  defined  by  r =10' 5 and  q = 10+6. 


Figure  7.8.  Total  input  energy  time-histories  for  the  five-story  buildings  with  fixed-base,  base- 
isolation,  supplemental  damping,  and  semi-active  control  subjected  to  the  SOOE  component  of  El 
Centro,  the  Imperial  Valley,  CA  earthquake  of  May  18,  1940  scaled  to  an  EPA  = 0.4g.  Parameters 
used  for  Bouc-Wen  Model  for  the  stories  are  A = 1,  y = 0.5,  n = 20,  a = 0.02,  and  for  the  base 
isolation  system  are  A = 1,  fi  = y = 0.5,  n = 3,  a = 0.6.  The  total  damping  of  the  building  with 
supplemental  damping  is  40%.  The  weighting  matrices  for  semi-active  building  are  defined  by  r =10* 
5 and  q = 10+6. 
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8.  CONCLUSIONS 


Using  20  accelerograms,  the  distribution  of  earthquake  input  energy  in  a one  and  a five- 
story  building  with  fixed  base,  based  isolation  system,  and  semi-active  control  system  are 
examined.  The  following  observations  and  conclusions  are  presented: 


8.1.  SDOF  structures: 

1.  Earlier  studies  on  a few  accelerograms  indicated  that  structural  properties: 
ductility,  damping,  and  the  post-  to  pre-yield  stiffness  ratio  for  bilinear  behavior 
do  not  have  a significant  influence  on  the  earthquake  input  energy.  Using  20 
accelerograms  this  study,  however,  shows  that  ductility,  damping,  and  the  shape 
of  the  hysteresis  loop  have  a major  influence  on  the  input  energy  spectra  and  its 
distribution  among  various  energy  components.  Results  indicate  that  in 
estimating  the  energy  parameters,  the  influences  of  the  structural  properties  are  as 
significant  as  the  influences  of  ground  motion  characteristics. 

2.  This  study  shows  that  as  the  effective  peak  acceleration  increases,  the  input 
energy  also  increases;  indicating  that  the  input  energy  is  related  to  the  intensity  of 
ground  motion.  The  energy  ratios  such  as  the  ratio  of  the  maximum  hysteretic 
energy  to  the  maximum  input  energy  EhJE,rm  are  not  affected  by  the  effective 
peak  acceleration.  Therefore,  scaling  an  accelerogram  does  not  change  the 
distribution  of  earthquake  input  energy  among  energy  components.  Results 
indicate  that  the  influence  of  the  duration  of  strong  motion  on  the  input  and 
hysteretic  energies  is  as  significant  as  the  influence  of  frequency  content, 
particularly  for  non-linear  structures.  As  the  ductility  ratio  increases,  the 
influence  of  the  duration  of  strong  motion  on  the  input  energy  spectra  becomes 
more  significant,  particularly  in  the  vicinity  of  the  predominant  period  Tpe. 

3.  For  the  20  accelerograms  used,  the  predominant  period  Tpa  (computed  from  the 
acceleration  spectra)  is  different  from  Tpe  (computed  from  the  energy  spectra) 
indicating  that  the  strength  and  energy  demands  do  not  necessarily  occur  at  the 
same  period.  Therefore,  in  developing  an  energy-based  design  approach  one 
should  examine  both  strength  and  energy  demands. 

4.  The  study  concludes  that  the  elastic  input  energy  spectra  E„  do  not  necessarily 
provide  a conservative  estimate  of  the  inelastic  input  energy  spectra. 

5.  Three  energy  ratios  based  on  hysteretic  energy  were  compared:  the  maximum 
ratio  of  the  hysteretic  to  input  energy  (Eh/Eir)m,  the  ratio  of  the  maximum 
hysteretic  energy  to  the  maximum  input  energy  Ehm/Ejrm , and  the  equivalent 
number  of  yield  excursions  Neq.  It  is  found  that  the  (Ef/Eir)m  ratio  reflects  the 
damage  potential  associated  with  the  largest  yield  excursion,  and  the  Ehm/Eirm  and 
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Neq  ratios  reflect  the  damage  potential  associated  with  the  total  number  of  yield 
excursions  and  the  cumulative  inelastic  deformation  for  the  entire  duration  of 
accelerogram.  The  study  shows  that  (Eh/Eir)m  is  independent  of  the  duration  of 
strong  motion  and  period  of  structure,  whereas  EhnJE,rm  is  independent  of  both 
only  for  periods  less  than  1 5.  Results  show  that  as  the  duration  of  strong  motion 
increases  the  equivalent  number  of  yield  excursions  Neq  increases;  indicating 
more  seismic  vulnerability.  Nonetheless,  one  cannot  make  a definitive 
conclusion  from  Neq  spectra  for  comparing  the  damage  potential  of  various 
structures. 

6.  It  is  observed  that  the  number  and  the  rate  of  yield  excursions  depend  on  the 
ground  motion  characteristics  (intensity,  frequency  content,  and  duration  of  strong 
motion),  and  structural  properties  (period  of  structure,  ductility,  damping,  and 
hysteretic  behavior).  The  damage  potential  for  the  largest  yield  excursion  is  as 
significant  as  the  total  yield  excursion  in  evaluating  the  damage  potential  of 
structures. 

7.  The  results  indicate  that  small  damping  ratios  (less  than  5%)  have  a minor 
influence  on  the  input  energy  and  a major  influence  on  the  damage  potential  of 
structures.  Therefore,  input  energy  should  not  be  considered  as  a sole  damage 
index.  For  damping  ratios  greater  than  5%,  however,  damping  reduces  the  input 
energy  in  the  vicinity  of  the  predominant  period  Tpe. 

8.  Earthquake  input  energy  reflects  the  intensity  of  earthquake  ground  motion  and  is 
more  appropriate  for  measuring  the  intensity  than  using  the  effective  peak 
acceleration  EPA  (equation  1.1),  Spectrum  Intensity  SI  (equation  1.2),  and 
Accelerogram  Intensity  Ie  (equation  4.6)  since  E„  considers  both  the 
characteristics  of  ground  motion  and  the  structural  properties.  The  maximum 
ordinate  of  the  elastic  input  energy  spectrum  with  5%  damping  may  be  used  for 
scaling  the  records  to  the  same  intensity. 

8.2.  Five-story  building  with  fixed  base,  base-isolation,  supplemental  damping, 
and  semi-active  control: 

1.  For  a given  effective  peak  acceleration,  the  distribution  of  energy  and  story  shear 
through  the  height  of  the  five-story  building  (with  fixed  base,  base-isolation, 
supplemental  damping,  and  semi-active  control)  is  mostly  independent  of  the 
duration  of  strong  motion  and  the  frequency  content  of  the  accelerograms. 

2.  Total  drift  and  total  input  energy  for  the  building  with  base-isolation, 
supplemental  damping,  and  semi-active  control  are  significantly  less  than  those 
for  the  fixed-base  building. 
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3.  The  ductility  and  hysteretic  energy  demands  do  not  necessarily  occur  in  the  same 
story.  The  story  shear  is  the  largest  in  the  first  story,  whereas  ductility  and  energy 
dissipation  may  not  be  the  largest  in  the  first  story. 

4.  Base-isolation  is  effective  in  reducing  the  damage  potential  of  the  5-story 
structure  by  dissipating  most  of  the  earthquake  input  energy  in  the  base  isolation 
system.  Nonetheless,  if  the  period  of  the  base-isolation  system  is  close  to  the 
predominant  period  of  the  ground  motion,  the  hysteretic  energy  demand  and 
consequently  the  damage  potential  for  the  building  with  base-isolation  will  be 
larger  than  those  for  the  building  with  fixed  base.  The  predominant  period  of  the 
ground  motion  primarily  depends  on  the  soil  conditions.  It  also  depends  on  the 
fault  and  path  characteristics  and  epicentral  distance  (specially  in  near  fault 
regions). 

5.  For  the  5-story  building  with  fixed  base  subjected  to  the  SOOE  component  of  the 
El  Centro  accelerogram,  as  the  yield  deformation  decreases,  base  shear  decreases 
significantly,  and  the  total  input  energy  remains  approximately  the  same.  It  is 
found  that  ductility  has  a pronounced  influence  on  the  distribution  of  input  energy 
among  energy  components. 
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APPENDIX  A.  BOUC-WEN  MODEL 


Different  models  are  used  to  characterize  the  nonlinear  behavior  of  structures  such  as 
Ramberg-Osgood  (Ramberg  and  Osgood,  1943),  Bouc-Wen  (Bouc,  1967;  Wen,  1976  and 
1980),  Penzien  and  Liu  (Penzien  and  Liu,  1969),  Takeda  (Takeda,  1970),  Sina  (Saiidi  and 
Sozen,  1979),  and  Q-Histeresis  (Saiidi  and  Sozen,  1979),  Otani  (Otani,  1981), 
Vinodgradov-Pivovarov  (Vinodgradov  and  Pivovarov,  1986).  The  Bouc-Wen  model 
which  is  extensively  used  by  researchers  (see  Wen,  1989)  to  represent  the  force- 
deformation  relationship  is  used  in  this  study. 

For  the  Bouc-Wen  model.  Figure  A.l,  the  equation  of  motion  can  be  written  as 


(A.l) 


where 

u = displacement, 
uy  = yield  displacement, 

iig  = ground  acceleration, 

£0=  viscous  damping  ratio, 

co0  = pre-yield  natural  frequency, 

a = post-  to  pre-yield  stiffness  ratio. 

In  the  above  equation,  z is  a non-dimensional  hysteretic  parameter  which  is  related  to 
displacement  u through  the  following  first-order  non-linear  differential  equation 


(A.2) 


Parameters  y and  /?  define  the  shape  of  the  hysteretic  loop  (softening  or  hardening); 
parameter  A controls  the  restoring  force  amplitude  and  tangent  stiffness,  and  n defines  the 
smoothness  of  the  transition  from  elastic  to  inelastic  range  in  the  force-deformation 
relationship.  This  model  can  be  used  for  a variety  of  materials.  For  example  A = 1,  y = 
P = 0.5,  n = 20  and  a = 0 can  approximately  represent  the  elasto-plastic  behavior. 
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Figure  A.l  Single  degree  of  freedom  system  used  for  Bouc-Wen  model 
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APPENDIX  B.  EARTHQUAKE  INPUT  ENERGY  FOR  MDOF 

STRUCTURES 


The  equation  of  motion  for  an  n degree  of  freedom  shear  building  (Figure  7.1. a)  can  be 
written  as 

Ft  (0  + Fd  (t)  + Fs  (0  = -Mr  ug  (t)  (B.  1 ) 

where  Ff  is  an  n -dimensional  inertia  force  vector,  Fd  is  an  n -dimensional  damping 

force  vector,  Fs  is  an  n -dimensional  non-linear  stiffness  force  vector,  M is  the  mass 
matrix,  r is  an  n -dimensional  influence  vector,  and  iig  denotes  ground  acceleration. 

Transposing  the  matrix  equation  (B.l),  multiplying  both  sides  by  du  = ii  dt  and 
integrating,  we  obtain  the  energy  balance  equation  as 

£>/  du  + JV/  (0  du  + JV/  (0  du  = - JVr  ug  (/)  du  (B.2) 

where  u is  an  n -dimensional  relative  displacement  vector.  The  relative  kinetic  energy  is 


Ekr  = J Mw(0  du  = V2  u1  M ii 

The  damping  energy  is  given  by 

Ed  = j'0F/(t)udt 
The  absorbed  energy  vector  is  given  by 

£„  = {//(  t)udt  = E„+Es 

the  elastic  strain  energy  for  the y'th  story  of  building  is 


(B.3) 


(B.4) 


(B.5) 


(B.6) 
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where  Q ■ and  k t are  the  shear  force  and  elastic  stiffness  of  the  y'th  story,  respectively. 

Once  the  elastic  strain  energy  is  obtained,  the  hysteretic  energy  can  be  computed  from 
equation  (B.5).  Finally,  the  relative  input  energy  is  obtained  as 


E 


ir 


Mr  iig  ( t ) ii  dt 


(B.7) 


Therefore,  the  vector  equation  (B.2)  can  be  written  as 


Ekr  + Ed  + E,  + Eh  - Eir 


(B.8) 
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APPENDIX  C.  LIST  OF  SYMBOLS 


a,j  A regression  parameter  depending  on  the  ductility  and  the  duration  of 

strong  motion 

A Parameter  of  Bouc-Wen  model  describing  tangent  stiffness 

bu  A regression  parameter  depending  on  the  ductility  and  the  duration  of 

strong  motion 

B Control  input  matrix 

c Viscous  damping  coefficient  of  SDOF  structure 

Cb  Viscous  damping  coefficient  for  base-isolation  system 

c,j.  A regression  parameter  depending  on  the  ductility  and  the  duration  of 
strong  motion 

C Viscous  damping  coefficient  of  each  story  of  the  structure 

d/j  A regression  parameter  depending  on  the  ductility  and  the  duration  of 

strong  motion 

Ea  Absorbed  energy 

Ec  Energy  corresponding  to  energy 

Ed  Equivalent  viscous  damping  energy 

Edm  Maximum  equivalent  viscous  damping  energy 

Edis  Energy  dissipated  through  damping  and  hysteretic  behavior 

Edism  Maximum  dissipated  energy  through  damping  and  hysteretic  behavior 

Edist  Total  dissipated  energy  through  damping  and  hysteretic  behavior 

Edt  Total  equivalent  viscous  damping  energy 

Eh  Hysteretic  (yielding  or  plastic  strain)  energy 

Ehm  Maximum  hysteretic  energy 

Eht  Total  hysteretic  energy 

Eia  Absolute  input  energy 

Eir  Relative  input  energy 

Eirm  Maximum  relative  input  energy 

Eir,  Total  relative  input  energy 

Ekr  Relative  kinetic  energy 

Ekrm  Maximum  relative  kinetic  energy 

Ekr,  Total  relative  kinetic  energy 

Es  Elastic  strain  energy 

Esm  Maximum  elastic  strain  energy 

Est  Total  elastic  strain  energy 

EPA  Effective  peak  ground  acceleration 

fd  Damping  force 

f.  Inertia  force 

fs  Restoring  force 

Fc  Control  force 

Fd  Damping  force 

Fj  Inertia  force 

Fs  Restoring  force 
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Fy 

FS 

g 

G 

H 

I 

h 

G 

J 

k 

h 

K 

m 

mb 

n 

Neq 

P 

PGA 

PGV 

q 

Q 

Q* 

r 

R 

Sv 

SI 

tdb 

tdi 


Th 

Tg 

T„ 


ll 

ug 

Uni 
Ur  el 
Uy 
Uyl 

ii 

iit 

VE 


Yield  force 
Fourier  Spectra 
Gravity  acceleration  of  ground 
Feedback  gain  matrix 
Disturbance  input  vector 
Unit  matrix 

Intensity  of  accelerogram 

Seismic  Index 

Performance  index 

.Stiffness  of  the  SDOF  structure 

.Stiffness  of  the  base  isolation  system 

.Stiffness  for  each  story  of  the  structure 

Mass  of  the  SDOF  structure 

Mass  of  the  base  for  a base-isolated  building 

A parameter  for  Bouc-Wen  model  controlling  the  smoothness  of  transition 

from  elastic  to  inelastic  region  of  the  force-displacement  curve 

Equivalent  number  of  yield  excursions 

A matrix  for  the  solution  of  Riccati  equation 

Peak  ground  acceleration 

Peak  ground  velocity 

A scalar 

Story  shear;  a weight  matrix  in  Chapter  7 
Base  shear 
A scalar 

Response  modification  factor;  a weight  matrix  in  Chapter  7 

Spectral  velocity 

Spectrum  Intensity 

Bracketed  duration  of  strong  motion 

Intensity-based  duration  of  strong  motion  defined  by  Trifunac  and  Brady 
(1975) 

Natural  period  of  the  base  isolation  system 
Predominant  period  of  the  earthquake  ground  motion 
Fundamental  period  of  the  structure 

Predominant  period  of  earthquake  ground  motion  based  on  acceleration 
response  spectra 

Predominant  period  of  earthquake  ground  motion  based  on  acceleration 
energy  spectra 

Relative  displacement;  deformation 

Ground  displacement 

Maximum  deformation 

Relative  displacement 

Yield  deformation 

Yield  deformation  for  the  first  story 

Relative  velocity  of  the  structure  with  respect  to  ground 

Absolute  (total)  acceleration  of  the  structure  with  respect  to  ground 

Equivalent  velocity 
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W Weight  of  the  each  floor  of  the  structure 

z A non-dimensional  hysteretic  parameter 

a Post-  to  pre-  yield  stiffness  ratio 

P Parameter  of  Bouc-Wen  model  defining  the  shape  of  hysteresis  loop 

y Parameter  of  Bouc-Wen  model  defining  the  shape  of  hysteresis  loop 

fu  Ductility  ratio 

Mmax  Maximum  inter-story  ductility  ratio 

% Damping  ratio;  a mass  vector  in  Chapter  7 

% b Damping  ratio  for  base  isolation  system 

co0  Pre-yield  frequency  for  Bouc-Wen  model 

Viscous  damping  ratio  for  Bouc-Wen  model 
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